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Abstract 

A quantum-kinetic formulation of the dynamical evolution of a high-energy non-equilibrium gluon system at finite 
density is developed, to study the interplay between quantum fluctuations of high-momentum (hard) gluons and 
the low-momentum (soft) mean color-field that is induced by the collective motion of the hard particles. From the 
exact field-equations of motion of QCD, a self-consistent set of approximate quantum- kinetic equations are derived by 
separating hard and soft dynamics and choosing a convenient axial-type gauge. This set of master equations describes 
the momentum space evolution of the individual hard quanta, the space-time development of the ensemble of hard 
gluons, and the generation of the soft mean-field by the current of the hard particles. The quantum-kinetic equations 
are approximately solved to order g 2 (l + gA) for a specific example, namely the scenario of a high-energy gluon beam 
along the lightcone, demonstrating the practical applicability of the approach. 



I. INTRODUCTION AND SUMMARY 



The physics of high-density QCD becomes an increasingly popular object of research, both from the experimental, 
phcnomenological interest, and from the theoretical, fundamental point of view. Presently, and in the near future, 
the collider facilities HERA (ep, eAl), Tevatron (pp, pA), RHIC and LHC (pp, AA) are able to probe new regimes 
of dense quark-gluon matter at very small Bjorken-x or/and at large A, with rather different dynamical properties. 
The common feature of high-density QCD matter that can be produced in these experiments, is an expected novel 
exhibition of the interplay between the high-momentum (short-distance) perturbative regime and the low-momentum 
(long-wavelength) non-perturbative physics. For example, with HERA and Tevatron experiments, one hopes to gain 
insight into problems concerning the saturation of the strong rise of the proton structure functions at small Bjorken-z, 
possibly due to color-screening effects that are associated with the overlappping of a large number of small- a; partons. 
Another example is the anticipated formation of a quark-gluon plasma in RHIC and LHC heavy ion collisions, where 
multiple parton rescattering and cascading may generate a high-density environment, in which the collective motion 
of the quanta can give rise to non-abelian long-wavelength excitations and screening of color charges. 

In any case, the study of coherent low-momentum excitations in QCD, that are generated by, and interacting with, 
the high-momentum partonic color charges, is of fundamental interest in several respects: Firstly, it provides insight 
into the basic features of non-abelian multiparticle dynamics and a step towards a rigorous decription of parton trans- 
port properties in a dense environment. Secondly, it may help to resolve current problems encountered in perturbative 
QCD, for instance the absence of static magnetic color-screening Q, the problem of infrared renormalons || connected 
with the resummation of perturbation theory in the small- a; regime, or, the problem of confinement associated with 
collective 'glue'-behavior of non-perturbative gluons ||. Interesting progress in these areas is continously being made, 
and consistent schemes have emerged to perform calculations of the parton evolution at very small- x |Q] , at very large 
density |5||| and for high-temperature QCD of a quark-gluon plasma 0. 

Most progress in the context of bulk multi-parton dynamics at high density has been made by studying 'hot 
QCD' with a thermally equilibrated quark-gluon system at very high temperature T. 'Hot QCD' has the attractive 
advantage that the parton density is homogenous and isotropic in momentum, and its exact form oc T 3 is known, 
since T»Asi 200 MeV is the only energy scale in the problem. For this academic scenario, inconsistencies of former 
perturbative calculations have been resolved by gauge-invariant resummation techniques || as studied in various 
applications |J, and moreover, a self-consistent kinetic theory has been formulated |Tc| |. 

The present paper, extending previous work of Ref. jtl), is to be viewed in this very context: it takes the 'hot 
QCD' developments as inspirational guideline, but aims to describe the opposite physics extreme, namely a highly 
non-equilibrium Q non-uniform and non-isotropic parton system. Specifically, the attempt is made to derive from 
first principles a self-consistent kinetic description for a non-quilibrium scenario of a gluon beam directed along the 
lightcone, that is, a high-density system of gluons, moving with very large energies ka ~ k z 3> k± 3> A along a beam 
direction (the fc z -axis), as it would be typical for the initial stage of a high-energy collider experiment (an extreme 
example is a collision of two heavy nuclei at the LHC, involving many thousands of gluons coming down the beam 
pipe). For simplicity the quark degrees of freedom are ignored, but are straightforward to include. 

As illustrated schematically in Fig. 1, the initial multi-gluon state is imagined as a highly Lorentz contracted sheet 
of bare gluons, characterized by a very large momentum scale Q (e.g. in an ultra-relativistic nuclear collision, the 
typical momentum transfer of hard scatterings that materialize the gluons out of the colliding beam nuclei) . Hence the 
typical energy and longitudinal momentum of the initial gluons is ~ Q. The subsequent evolution of these bare quanta 
is, at leading order a s , well known to lead to a rapid multiplication and diffusion of gluons through real and virtual 
radiation, corresponding to brems-strahlung and Coulomb-field regeneration, respectively |T^ |. As a consequence, 
the typical gluon momenta both in longitudinal and transverse direction, decrease (see Fig. la). As long as the 
average transverse momentum is sufficiently large, k± > pt ~ 1 — 2 GeV, a s (fi ) < 1, a perturbative description of 
the evolution of the gluon density G is appropriate, but when k± < /i, non-perturbative dynamics is expected to take 
over, governed by the collective infrared behavior of a large number of long-wavelength gluons. If the number density 
of low-mometum gluons below /i is large, their dynamics may approximately be described classically [|[[l3) in terms 
of a coherent mean field A (see Fig. lb). 



* The term 'non-equilibrium' is used in the sense of statistical many-body physics, describing a quantum system far off the 
state of maximum entropy and thermal equilibrium. Such a non-equilibrium system may in general be spatially inhomogenous 
and anisotropic in momentum, in contrast to a homogenous, thermal ensemble, or translation invariant system in vacuum. 



1 





FIG. 1. Non-equilibrium scenario of gluon beam along the lightcone: a) The initial multi-gluon state, prepared at time 
to — at the hard scale Q with initial condition Ao(0, Q + , Q 2 ), develops forward in time which is described by the evolution 
of the gluon propagator A(r, K + , K\) being a function of both space-time r M = (t,r) and momentum — (E,K). The 
gluons, propagating with large K + — E + K z S> K± along the z-axis, are accompanied by real and virtual radiation which 
causes a diffusion in both transverse direction r± and transverse momentum K± as time goes on: The emission of gluons 
increases the multiplicity and decreases the average transverse momenta (Kj_) at given lightcone time r~ = t — z and position 
lightcone position r + = t + z. b) Top: Qualitative picture of time evolution of the typical transverse momentum {K\) of hard 
gluons, where the earliest emitted daughter gluons have the largest K\ < Q 2 and later produced gluons have much smaller 
K\ . Eventually modes with K\ < fi 2 will be populated significantly. Bottom: Corresponding time development of the number 
density G of hard gluons from initial value Go and of the average soft field A that is induced by the population of gluons with 
K\ < fi 2 , starting from zero initial value. Speculatively, one would expect a saturation at asymptotic times due to screening 
of further gluon emission by the presence of the soft mean field. 
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Given this heuristic picture, the near-at-hand rationale is therefore to subdivide the dynamical development of 
the gluon ensemble into a perturbative quantum evolution in the short-distance regime Q 2 > fei > M 2 , and a non- 
perturbative, but classical, mean-field in the long- wavelength regime k\ < p? . The corresponding degrees of freedom 
are referred to as hard gluons for fcj_ > /i, whereas excitations with k± < p, represent the soft mean field. 

Because the hard gluons have small transverse extent A ~ l/fc_i_ < 0.2 fm (for p, = lGeV), they can be considered, 
locally in space-time, as incoherent self-interacting quanta, if the interparticle distance is significantly larger than 
A. On the other hand, when the typical transverse momenta drop below p, the gluons begin to act coherently, 
and collectivity arises, because the motion taking place over a distance scale 1//* or larger, involves coherently a 
large number of hard particles, which gives rise to an average soft color field. The crucial point of this hard-soft 
separation is that the over long distances A > l/p, the soft mean field represents the average gluon motion, but at 
short-distances A < 1/p, the hard gluons may be described approximately as in free space. Certainly, such a rigid 
division of hard and soft physics in terms of a single parameter p, is at his point an arbitrary and idealizing definition. 
However, the arbitrariness can in principle be removed by considering the variation with respect to p, as in the 
usual renormalization-group framework. This interesting task is beyond the scope of this paper, and remains to be 
addressed in the future. 

The non-quilibrium scenario of a lightcone beam of gluons along the lightcone has two major advantages over 
the opposite thermal equilibrium extreme, the isotropic quark-gluon plasma: First, it favors the 2-scale separation 
between hard and soft physics. Second, it allows to choose an axial-type gauge which eliminates to large extent the 
problems of non-linearities and of ghost degrees of freedom that are encountered in usual covariant gauges. The 
2-scale separation arises naturally here, because Lorentz contraction and time dilation along the beam direction plus 
the limited transverse momenta, force the hard gluon fluctuations and self-interactions to be highly localized to short- 
distances, and separates their quantum motion from the low-momentum mean-field dynamics over comparably long 
distances. On the other hand, the choice of a non-covariant axial gauge, characterized by a directed four-vector n 
along a fixed axis, is very suggestive, because the geometry and kinematics allows to choose n paralell to the gluon 
momentum k z , in which case perturbative QCD calculations formally reduce in many respects to the abelian QED 
counter parts. This is not possible for an isotropic thermal system, where all possible directions of gluon motion 
are equally probable. Given these premises, the quantum dynamics is dominated by the self-interactions of the hard 
gluons, which make them fluctuate localized around the lightcone, whereas the kinetic dynamics can well be desribed 
statistical-mechanically in terms of mutual interactions among them and in the presence of their generated soft mean 
field. As elaborated in Ref. pH , these notions are the keys to formulating a quantum-kinetic description, by combining 
standard techniques of parton evolution and renormalization group, with relativistic many-body transport theory 

The main result of this study within the outlined physics framework, is a set of three master equations, which couple 
the quantum evolution of short-distance fluctuations of the individual hard gluons, the space-time development of the 
gluon system as-a-whole, and the generation of the soft mean field: 

(i) an evolution equation for the spectral density p of each individual hard gluon, which determines the intrinsic gluon 

distribution of a hard gluon in accord with mass- and coupling-constant renormalization, and which dresses up 
the bare initial gluons to renormalized 'quasi-particles'. 

(ii) a transport equation for the space-time development of the whole ensemble of these renormalized gluons with 

respect to their propagation in the self-generated soft mean field, as well as due to their scatterings off each 
other, which determines the physical gluon phase-space density G. 

(iii) a Yang-Mills equation for the generation of the soft mean field A, which is induced by the effective color current 
of the hard, renormalized gluons, where the current is obtained from the momentum-weighted gluon phase-space 
density. 

Although this set of equations appears at first sight to be of impractical complexity, it allows in fact for a practical 
applicable calculation scheme, as will be demonstrated with an explicit sample calculation. 

To arrive at the above master equations, three essential aspects of the problem have to be merged: first, the physics- 
dictated aspect of space-time, kinematics and geometry, second, the quantum field aspect of gluon excitations and 
self-interactions, and third, the statistical aspect of multi-particle interactions in the presence of the mean field. The 
non-trivial interconnection of these aspects require to work directly at the level of equations of motion, rather than on 
the level of Feynman diagrams, because the relative proportions and interactions of hard and soft quanta must can 
only be calculated self-consistently from the equations of motion. 

The strategy for deriving the above three-some of master equations follows closely the previous work of Refs. JlT| . 
The path-integral representation of the Yang-Mills action gives an infinite set of equations of motion for the non- 
equilibrium n-point Green functions, which is the well known analogue of the BBGKY hierarchy This hierarchy, 
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which represents the exact theory, is truncated to a system of equations involving only the 1- and 2-point functions, by 
arguing that higher-order correlators n > 3 are comparably small. To achieve self-consistency of the truncated set of 
equations at the n = 2 level, the n > 3 functions must be implicitely lumped into the 1- and 2-point functions. After 
separating hard and soft field modes, as alluded before, the 1-point function is identified with the soft average field 
Afj_ = (A^(x)) and the 2-point function is given by the hard gluon correlator iA^ — (a fl (x) a v (y))p, where and 
represent the soft and hard modes, repectively. The truncated set of equations of motion then involves the non- 
equilibrium version of the Dyson-Schwinger equation for A and the classical Yang-Mills equation for the soft mean- field 
A. The two field-equations of motion for A and A can be cast into much simpler quantum-kinetic equations with 
the help of the Wigner-function technique and gradient expansion, and the assumption of 2-scale separation implying 
that the long-wavelength A-ficld is slowly varying on the short-distance scale of the hard quantum fluctuations. The 
result is then the above set of master equations. 



A powerful theoretical framework to derive from the exact field equations of motions the above approximate quantum 
kinetic equations, is the so-called Closed-Time-Path formalism (CTP). The CTP formalism is a general tool for 
treating initial value problems of irreversible multi-particle dynamics in quantum field theory. It therefore provides an 
appropriate language to describe the problem of non-equilibrium gluon dynamics within a well-established theoretical 



framework. Originally introduced by Schwinger Jlf 
is documented in great detail in the literature flTi 



and Keldysh [|16| the CTP formalism and its diverse applications 
In particular, I refer to Ref. UM, where the CTP method is 



applied to high-energy QCD, and to Appendices B and C. 
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FIG. 2. Difference between in-out formalism of usual quantum field theory (in free space or 'vacuum') and the in-in 
formalism of the CTP formulation (in the presence of surrounding particles or 'medium'), a) The in-out amplitude describes by 
the evolution of an asymptotic \in) state at to — * — oo to an asymptotic \out) state at too — * oo by means of the time evolution 
operator V (to, too). Because U (to, too) = £^(£00,^0), forward and backward evolution are identical, and there is no correlation 
between the two time branches. Consequently, the Feynman propagator Gf = contains the full dynamics of the 2-point 
correlations, b) The in-in amplitude starts at to with a non-trivial initial multi-particle state described by the density matrix 
f>(to) and evolves again by means of time evolution operator U(to, too)- Due to statistical interactions among the many evolving 
particles acting as a medium, in this case [>(to)U(to,too) 7^ (too, to)- Consequently, the statistical correlation between the 
two time branches have the effect that Gf 7^ G^, and moreover require the introduction of additional correlation functions G> 
and G< to account for the cross-talk between upper and lower time branches. 
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The fundamental starting point of non-equilibrium field theory in the CTP formalism is to write down the in-in 
amplitude Zp for the evolution of the initial quantum state \in) forward in time into the remote future. As reviewed 
in Appendix B, this generalizes the usual quantum field theory approach based on the vacuum- vacuum transition 
amplitude, or in- out amplitude, to account for the a-priori-prcscnce of medium particles described by the density 
matrix p(to) and to evolve this non-trivial initial state in the presence of the medium from to to t^ in the future, 
The in-in amplitude Zp is graphically depicted in Fig. 2, and formally it is given by Zp\J,p\ — ( in \ in )j 
where J = (^J + ,^J~) is an external source with components on the upper + and lower — time branch, and p(to) 
denotes the the initial state density matrix. From the path-integral representation of Zp one obtains then the non- 
equilibrium Green functions. The convenient feature of this Green function formalism on the closed-time path P is 
that it is formally analogous to standard quantum field theory, based on the vacuum- vacuum, or in- out amplitude 
Z[J] = ( out | in) j = (0\0)j, except for the fact that in the CTP formalism, the fields have contributions from 
both time branches. For details I refer to Appendix B, where the basics of the CTP formalism are summarized, in 
particular, how to obtain the path-integral for Zp that generates the Green functions on the closed-time-path P. 

The interpretation of this formal apparatus for the evolution along the closed-time path P is rather simple: If the 
initial state is the vacuum itself, that is, the absence of a medium generated by other particles, then the density matrix 
p is diagonal and one has \in) — > |0). In this case the evolution along the + branch is identical to the anti-time ordered 
evolution along the — branch (modulo an irrelevant phase) , and space-time points on different branches cannot cross- 
talk. In the presence of a medium however, the density matrix contains off-diagonal elements, and there are statistical 
correlations between the quantum system and the medium particles (e.g. scatterings) that lead to correlations between 
space-time points on the + branch with space-time points on the — branch. Hence, when addressing the evolution 
of a multi-particle system, both the deterministic self-interaction of the quanta, i.e. the time- (anti-time-) ordered 
evolution along the + (— ) branch, and the statistical mutual interaction with each other, i.e. the non-time-ordered 
cross-talk between + and — branch, must be included in a self-constistent manner. The CTP method achieves this 
through the time integtation along the contour P. Although for physical observables the time values are on the + 
branch, both + and — branches will come into play at intermediate steps in a self-consistent calculation. 

The outline of the paper is as follows: In Section 2 the field equations of motion for the hard gluon propagator 
and the soft mean field are derived from the path-integral representation of the in-in amplitude Zp for non-covariant 
gauges. After separating hard and soft degrees of freedom, two key approximations are made, that allow to cast the 
infinite hierarchy of exact equations of motion in terms of a truncated system of only two approximate equations, 
namely a Dyson-Schwinger equation for the hard gluon propagator, and a Yang-Mills equation for the soft field. In 
Section 3 the transition to a quantum kinetic description is worked out. This requires one further key approximation 
in conjunction with a clear definition of quantum and kinetic space-time regimes, such that the aforementioned 2- 
scale separation is guaranteed. This also defines the limits for the applicability of the quantum kinetic approximation. 
Provided that the separability condition is satisfied, one finally arrives at the set of master equations discussed above, 
for which a systematic calculation scheme is proposed. In Section 4 an explicit calculation to solve the master equations 
is presented for the physics scenario depicted in Fig. 1. I consider the evolution of an initial incoherent ensemble 
of bare gluons moving colinearly along the lightcone as it proceeds in its momentum- and space-time development 
and generates its soft mean field. To avoid overkill of too many technical details, each Section is accompanied by 
Appendices. Appendix A defines the notation and conventions used throughout the paper. Appendix B reviews the 
basics of the CTP formalism. Appendix C discusses the application of the CTP method to QCD for non-covariant 
gauges. Appendix D shows the advantagous absence of ghosts in non-covariant gauges. Appendix E gives details on 
how to obtain from the in-in amplitude Zp an approximate effective action functional from which the of motion for 
the hard gluon propagator and the soft field are derived. Appendix F summarizes some basic analycity properties of 
the free-field propagators in the CTP formalism and discusses their relation to the gluon phase-space density. 

Finally some remarks on most closely related work in the literature (for an extended discussion, see introduction of 
Ref. @): 

Blaizot and Iancu pi have in a scries of papers developed a kinetic theory for 'hot QCD', i.e. the case of a high- 
temperature quark-gluon plasma. One of the key elements of their approach is the formulation of well-defined and 
consistent approximation scheme. I adopted many features of this approach to the present, rather different physics 
context. New is here the inclusion of the aspect of quantum evolution and renormalization. 

McLerran, Venugopalan, et al. [p4| , as well as Makhlin [ p5[ , have developed different approaches to calculate the 
quantum evolution of parton systems with lightcone dominance, i.e. in a beam-type scenario as considered in the 
present work. The McLerran- Venugopalan model also gives a predictive estimate for the feedback effect of the coherent 
mean field on the hard gluon evolution that generates this field. In several respects I follow a similar route. New in 
the present work is that it embodies in addition the aspect of space-time development of the evolution. 



5 



Boyanovski et al. have intensly studied the non-equilibrium evolution in scalar field theory, using extensively the 
techniques of the CTP formalism in conjunction with a large- TV expansion. Although the focus on this paper is rather 
different, many of the concepts and results in their papers concerning the first-principle time evolution of quantum 
system with associated particle production, dissipation, mean-field dynamics, etc., may serve as a scalar toy model 
for QCD. 



II. INTERPLAY OF 'HARD' AND 'SOFT' GLUON DYNAMICS 



A. The in-in amplitude for QCD in non-covariant gauges and the concept of approximation 

The in-in amplitude Zp introduced in Sec. 1 admits a path-integral representation which is the generating functional 
for the non-equilibrium Green functions defined on closed-time-path P, as discussed in Appendices B and C: 

Z P [K] = J VA expji (l[A,K,fj} , (1) 

where A® = (A^ + ,A^~) has two components, living on the upper (+) and lower (— ) time branches of Fig. 2, with 
DA = Y[a a -^^u + P*4^ ~ , and JC represents the presence of external sources. I consider here the class of non-covariant 



gauges defined by 27 2 



(n»A%{x)) =0, (2) 

where n M is a constant 4- vector, being either space-like (n 2 < 0), time-like (n 2 > 0), or light-like (n 2 = 0). The 
particular choice of the vector n M is usually dictated by the physics or computational convenience, and distinguishes 
axial gauge (n 2 < 0), temporal gauge (n 2 > 0), and lightcone gauge (n 2 = 0). Referring to Appendix D, the great 
advantage of these gauges is that the Fadeev-Popov ghosts decouple, so that in practical calculations the ghost degrees 
of freedom can be ignored, just as in abelian gauge theories. 

Then the action / in the exponential of (|l|) is given by (c.f. Appendix C), 

/ [A,K] = Iym [A] + I G f [n-A] + JC[A] , (3) 

containing the Yang-Mills action Iym, the gauge fixing term Igf, and the initial state source term /C, containing 
multi-point correlations concentrated at t = to: 

Iym [A] = -\ [ d A xF«{x)F^ a {x) , Igf [n ■ A] = ~ [ d 4 x [n ■ A a (x)} 2 , (4) 

4jp 2a Jp 

JC[A]=JC {0) + [ d A xJC% a {x) A^ a (x) + \f d i xd i ylC < -%t b {x,y)A^ a (x)A^ b {y) + ... . (5) 
Jp 2 Jp 

The exact knowledge of the in-in amplitude Zp from ([!]), would require to the calculation of all Green functions 
up to infinite order, and would correspond to the full solution of QCD in non-equilibrium media. Rather than that, 
the realistic goal is to formulate a practical calculation scheme for the kinetic evolution of a multi-guon system. In 
order to make progress, one needs to make reasonable approximations that are consistent with the specific physical 
problem under study, and truncate the infinite hierarchy of Green functions. 

In this Section a closed set of approximate equations is derived that are in principle solvable, given a suitable physics 
scenario. The basic idea is to describe an evolving gluon system in terms of two distinct components, namely, hard, 
short-range quantum fluctuations and soft, long-wavelength collective excitations, which I assume to be separable by 
a characteristic space-time distance. It is clear that the relative proportions and interactions of hard and soft degrees 
of freedom must be calculated self-consistently from the equations of motion. 

Starting from the in-in amplitude (0), the strategy of procedure is the following: 

1. The exact expression of the in-in amplitude Zp = exp(iM / p) is rewritten in terms of soft and hard field modes by 
splitting the gauge field A^ — + a^. Therefrom one obtains an infnite set of coupled equations for the Green 
functions. In order to reduce this to a finite system, I make 
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• approximation 1: The functional Wp = — ilnZp is expressed in terms of connected 1- and 2-point 
functions Q^ 1 ', alone by eliminating Q^ 71 ' for n > 3 as dynamical variables. Then the expectation 
values of and describe the induced soft mean field A^ and the hard (soft) correlation functions 

(IV). 

From the truncated functional Wp the corresponding effective action Tp is obtained, which generates the desired 
self-consistent equations of motion for A^, and D^. Here I make 

• approximation 2: It is assumed that the soft field dynamics can be treated classically by the non- 
propagating average field A^, and that the long-range propagation of soft modes, described by may 



be ignored at this level, i.e. Z? M „ <C A^A V . This assumption is motivated by the widely studied J24 13 
observation that a classical treatment of the long-distance dynamics of bosonic quantum fields at high 
density, obeying the classical field equations, should provide a good approximation, if the soft modes are 
sufficiently occupied. 

The original infinite equation system can then be reduced to a Yang-Mills equation for the classical, soft field A^, 
as it is induced by the current of hard quanta, and a Dyson-Schwinger equation for the hard propagator A„y subject 
to the presence of the soft mean field and to quantum fluctuations. These field equations of motion are still of very 
intractable non-linear character. They are further simplified to quantum-kinetic equations in Sec. III. 



B. Separating soft and hard dynamics 

The first step in the strategy is the separation of soft and hard physics in the path-integral formalism with Green 
functions of both the soft and hard quanta in the presence of the soft classical field that is induced by and feeding 
back to the quantum dynamics. A frequently used method for separate treatment of quantum and classical dynamics 
in field theory is the so-called 'background field method' |^9| which has been studied, e.g., in the context of dynamical 
symmetry breaking, vacuum structure, confinement and gravity, or for hot plasmas in finite temperature QCD. Within 
the background field method, one would split up the gauge field appearing in the classical action into an external 
classical background field and a quantum field which remains the sole dynamical variable in the path integral. I will 
however not follow this path, and rather prefer to treat soft and hard physics on equal footing, that is, to separate the 
gauge field into a soft classical field plus its soft quantum excitations, and a hard quantum field. Then both soft and 
hard fields can quantized and remain as dynamical variables a priori. 

The gauge field appearing in the classical action Iym [A] is split up into a soft ( long-range) part A^, and a 
hard (short-range) quantum field a^: 

A > ] = / (^ e+ * x -^ (fc) ^ 2 - fc ') + / -0ji e+ik ' XA U k ) °( k l-» 2 ) = W + (6) 

This is the formal definition of the terms 'soft' and 'hard', as used in this paper. The soft and hard physics are 
separated by the momentum scale /i which is at this point arbitrary. However, this arbitrariness can in principle 
be overcome by considering /j(x) as a dynamical variable depending on the space-time point x, rather than a fixed 
parameter, and determining it self-consistently from the local stability condition dA v {x) / dp? {x) = 0. From ([|) it is 
obvious, that the corresponding scale in space-time, A = l//i, divides soft and hard regimes in terms of the transverse 
wavelength of field modes, so that one may associate the soft field A^ being responsible for long range color collective 
effects, and the hard field embodying the short-range quantum dynamics. Consequently, the field strength tensor 
receives a soft part, a hard part, and a mixed contribution, 

= (f^[A\ + f«M + 4^[A,a]) (*) ■ (7) 



When quantizing this decomposed theory by writing down the appropriate m-m-amplitude Zp, one must be con- 
sistent with the gauge field decomposition (^) into soft and hard components and with the classical character of the 
former. Substituting the soft-hard mode decomposition (^) into ([!]), the functional integral of the in-in amplitude 



(Cll) becomes: 

Z P [K] = J VAVa exp{i(l[A] + I[a] + I [A, afj j , (8) 
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with the soft, hard, and mixed contribution, respectively, 



I [A] 



I [a] 



I \A, a] 



d A x 



d 4 x 



d 4 x 



4 M" 



J d A x K {1 l a A^ a + J d 4 xd 4 y A* a K, {2 \* b A"' ' 



(n-a°) 2 



2a 



+ d 4 xK, {1 }, a a^ a + / d i xd A y a^ a JC {2 },f a^ b + 



fJLV 



(9) 



(10) 



<t>i v r> a + F^r> a ) 

d 4 x{ -gf abc {S%K) (a^ b a^ c + A^ b a^ c + a^ b A^ c ) 

+ {dial) (Af 1 ' b A v ' c + a^ b A u > c + A»< b a v ' c ) 



^2 jace jbde 



ZA^AlaP^d"'* + A° il A b v A^ c a v - a + aja^a"' C A"' d 



(11) 



Note that in ( [U)| ) and (|TT|) terms involving 2-products oc a M A v do not contribute to -Zp, because their expectation 
value vanishes due to the soft-hard separation (|6[) which defines and A„ as complimentary. 

At this point I make approximation 1 from above. It is assumed that initial state can be represented as an ensemble 
of incoherent hard gluons, each of which has very small spatial extent Arj_ <A=l/ft corresponding to transverse 
momenta k\ ^> /j 2 . By definition of /z, the short-range character of these quantum fluctuations implies that the 
expectation value (a^) vanishes at all times. However, the long-range correlations of the eventually populated soft 
modes with small momenta k\ < /i 2 may lead to a collective mean field with non- vanishing (A M ). Accordingly, I 
impose the following condition on the expectation values of the fields: 



= for t < t 
> for t > t 



{ a a Jx) ) = for alU 



(12) 



Now I make approximation 2, that is, the quantum fluctuations of the soft field are ignored, assuming any multi-point 
correlations of soft fields to be small, 



( A£ (xt) ... AH (x n ) ) « ( A£ On) )...(A% (x n ) ) for all n > 2 , 

i.e. take A^ as a non-propagating and non-fluctuating, classical field. In particular, 

Wf v {x,y) = ( Al(x)Al(y) ) « < A^x) ) < A b u (y) ) , 

so that the limit " — ► 0" can be considered. 

As explained in more detail in Appendix E, the generating functional for the connected Green functions, 

W P [K] = -i In Z P [K] , 

which generates the infinite set of connected n-point Green functions via 
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i-i)G^:;^{xx,... iXn ) 



■W P [K] 



K=0 



( < ( Xl )... a% (x k ) A%+\ (x k+1 ) . . . A% (x n ) ) 



(c) 

p ' 



(13) 
(14) 

(15) 
(16) 



is truncated at level n > 3 on the basis of approximation (|12|) and (|14[). As a result, Wp becomes a functional of the 
1-point function (soft mean field A) and the 2-point function (hard propagator A) only: 



3» &%?(x,y) = {al{x)al{y))f = i&*(x,y) 



(17) 



These relations define the soft, classical mean field A, and the hard quantum propagator A in terms of expectation 
values of soft and hard field operators A M and a M , respectively. One now readily one obtains the effective action Tp 
(or proper vertex functional) via Legendre transformation (c.f. Appendix E): 
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Tp [Q] 



A, A 



= W F 



(18) 



which is a functional of only the soft field A and the hard propagator A as independent dynamical degrees of freedom. 

The equations of motion for the mean field A and for the hard propagator A in the presence of sources, follow now 
by differentiaition of ( |l8| ) with respect to A and A (c.f. Appendix E): 



ST 
STp 



Jf- = -/C«^(.x) - / d A yK ( - 2 ^ v ' ab (x,y)A u,b (y) 
„(x) Jp 



6A$(x,y) 



-K^^ ab (x,y) 



(19) 
(20) 



The self-consistent equations of motion of the dynamically evolving system are then obtained from ( |19[ ) and (|20j) by 
(i) imposing initial conditions in terms of the K, kernels at t — to, and (ii) by obtaining an explicit formula for Tp in 
terms of A^ and A M „. Concerning the initial conditions, I restrict to non-equilibrium initial states of Gaussian form 
(i.e. quadratic in the hard modes) and do not consider possible linear force terms. That is, I set 



IC w (x) 



0, IC^\x,y) 



> 



x°=y°=t 



(21) 



To obtain an explicit expression for Tp, the formal loop expansion of ( p"8] ) results in the well known Cornwall- Jackiw- 
Tomboulis formula iBOfl: 



A, A 



hff[A,a] - -Sp MA^A) 



+ r« 



A, A 



(22) 



where Sp[A£? . . .] = Tr f p d A X\d x-z ■ ■ ■ A{x\)B{x2) . . . stands for both the trace over color and Lorentz indices, as well 
as the integration over all space-time positions, hence giving the expectation value (AB . . .)p as defined in Appendix 
C. The physical interpretation of the various terms in this expression for Tp is the following pT|: 



(i) The first term is of order h and is given by the classical action (|])-(10) at A — A and switched-off sources JC: 



hff[A, a] 



I [A] + I [a] + I [A, a] 



A=A,K=0 



(23) 



Notice that in the limit a = 0, this reduces to the classical action for the soft mean field, I e ff[A, 0] = Iym [A] + 
I G F[n ■ A]. 

(ii) The second term in (|22|) is of order h 1 and contains the the contributions of the coupling between the 
soft mean field A and the hard quantum propagator A. The free propagator (see Fig. 3a) is given by 
[S 2 I e ff[A, a]/Sa(x)Sa(y)]A=Q-a=o with A switched off, which yields 



{A^XA^y) 



8 ab 8p(x,y) d^{d x ) d 2 , 



(24) 



where it is understood that the space-time arguments x and y in Ao satisfy (x — y)\ < 1//J, 2 , and 
d^{d x ) ee <T - n " d l + f d * + (n 2 + a-^) d * d * 



n ■ d x 



(n ■ d x y 



(25) 
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a ) ' 'free ' ' propagator: 

i A„ (x, v) = • • 

o y yJ/ x y 



b) "mean field" propagator: 



A 



i A (x, y) = 



x 



= i A (x, y) + 



A A 



A A 



c) "full" propagator: 

i A (x, y) = 



x 



= i A (x, y) + 



■0- + 



A<S> A<8> <8>A A<8> 

+ 6 + o 



d) ' full ' ' 3-point function: 



-ig V (x, y, z) = 





e) "full" 4-point function: 



-g W(x,y,z,w) = 




z w 



- X + v 




/) induced mean field: 
A (x) = 



= o + 



FIG. 3. Diagrammatics of the various terms used for the n-point functions appearing in the text. The 2-point function 
Cj( 2 ) = iA is the hard gluon propagator with the free-field propagator A (no interactions), the mean-field propagator A 
(including the interactions with the classical soft field A), and the full propagator A (including both mean-field and quantum 
(loop) interactions). Similarly, the connected 3-point function = —igV and the 4-point function t/' 4 ' = —g 2 W contain soft 
mean-field plus hard quantum contributions with internal full propagator A. Finally, the 1-point function is the soft mean field 
A that is generated by the hard gluons through the coupling to the full 3-point and 4-point functions V and W . 
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Even in the absence of quantum fluctuations, these contributions amount to a modification of the free propagator, 
such that the free propagator Ao becomes an effective propagator A in the mean field, dressed up by the presence 
of A. This mean field propagator (see Fig. 3b), denoted by A, is obtained from [S 2 I e ff[A, a] / 5a(x)Sa(y)] A=A -. a=0 

with finite A ^ 0, which results in 

(A-\l(x,y) = (V)#(*>V) -<^,J/), (26) 

where II denotes the self-energy contribution associated with the presence of the mean-field A ^ 0. Its explicit 
expression is given below in (|4q) . In other words, the effect of the mean field is to shift the pole in the free 
propagator Ao of ( |24| ) by a dynamically induced 'mass' term oc II, which can produce screening and damping 

effects. Note that Ap 1 = A It is important to realize that this mean field effect is still on the classical 

tree-level, and does not involve quantum fluctuations associated with radiative self-interactions among the hard 
gluons. 

(iii) The last term Tp^ in (22) represents the sum of all two-particle irreducible graphs of order h 2 , h 3 , . . . p0| , with 
the full propagator A, dressed by both the soft mean field and the quantum self-interactions (see Fig. 3c), 

Af u (x,y) ee Sfajx.y) + S&fajw) , (27) 

where the dependence of the full propagator on the soft mean field A is indicated by an explicit subscript, and 

. (28) 



^m^u - [a] 



A=0 <- A > V-v 



A=0 



Note that ^ A , that is, A^ denotes the full propagator for A — 0, whereas A is the free propagator 

(p4|). The real (dispersive) part of contains the virtual loop corrections associated with the gluon self- 
interactions, whereas the imaginary (dissipative) part contains the emission, absorption, and scattering processes 

(2) 

of hard gluons. In other words, T p embodies all the interesting quantum dynamics that is connected with 

( 2 \ 

renormalization group, entropy generation, dissipation, etc.. The explicit form of T p is diagramatically shown 
in Fig. 4, with the vertices and lines defined by Fig. 3. Suppressing color and Lorentz indices and employing a 
condensed notation, e.g., A(xi,X2) = A°'(ii, 12), the corresponding formula is, 



r (2) 
1 P 



A, A 



= r (1) + r (2) + r (3) + r (4) , (29) 



with the following contributions 
1 



r (1) =-^ J d^xd^y J d*x l d*y 1 W (x,y,xi,yi)A(yi,x 1 )A(y,x) (30) 

2 

r (2) = j^9 2 J d 4 xd 4 y J Jld^Xid^i V (x,x 1 ,x 2 ) A(xi,yi) A(x 2 ,yi) V{y 2 ,yi,y) A(y,x) (31) 

r( 3 ) = ^.9 4 J d A xd A y J '^fxid^yi Wo(x, xi, x 2 , x 3 ) A(x%, y%) A(x 2 , y 2 ) A(x 3 , ys) 

x W{y 3 ,y 2 ,y u y) A(y,x) (32) 

2 

T( 4 ) = J d 4 xd 4 y J Y[d 4 x i d 4 y i d' i Zi Wo(x,xi,x 2 ,xs) A{x 2 , z 2 ) A(a; 3 , z 3 ) 

xV{z 3 ,z 2 ,z 1 ) A(z 1 ,y 1 ) A(x 1 ,y 2 ) V{yi,y 2 ,y) A(y,x) . (33) 

The functions V and W are the full proper vertex functions for the 3-gluon and 4-gluon coupling, respectively. 
Their diagramatic representation is shown in Fig. 3d and 3e, and formally they are given by the functional 
derivatives of Tp at A ^ 0, namely, [5 n T p / Sa{x\) . . . 5a{xi)5 A(xi+\) . . . SA(x n )] A= -j. a=0 for n = 3 and n = 4, 
respectively: 
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-igV&(x,y,z) = -igV 0X %(x,y,z) + 0(g 3 ) 
-g 2 W?^{x,y,z,w) = ~g 2 W a x b ^(x,y,z,w) + 0{g 4 ) , (34) 

which, to lowest order in the coupling constant, reduce to the bare 3- and 4-gluon vertices Vo and Wo, respectively: 

Vo \%&V,z) = f abc {g^(dy-d x )uS i P {x 1 z)5 4 P (y,z) + g, w {d z - d y ) x 8%(y, x)S%(z, x) 

+ 



W 



gu\(d x - d z )^Sp(x,y)Sp(z,y)j (35) 

f;t(x,y,Z, W ) = - { ( <r oe f bde_ f ade f c b e^ + ( fabef cde _ fade fbcej 

+ (r e f dhe - r he r de ) gx^} 5%{ x , y )s 4 P ( z , w )s 4 P ( y , Z ) . (36) 



(2), 



FIG. 4. The 2-loop contribution , eqs 







, to the effective action Tp of eq. 



in the diagrammatic represen- 



ts 

tation of Fig. 3. Formally, F p is the sum of all two-particle irreduclible graphs with internal lines representing the full gluon 
propagators A and full 3- and 4-gluon vertices V and W. 



C. Equations of motion 

As scetched above and discussed in more detail in Appendix E, the equations of motion ( |l9| ) and (^) result from 
approximating the exact theory by truncation of the infinite hierarchy of equations for the n-point Green functions 
to the 1-point function (the soft mean field A(x)) and the 2-point function (the hard propagator iA(x,y)), with all 
higher-point functions being combinations of these and connected by the 3-gluon and 4-gluon vertices ~igV{x 1 y, z) 
and -~g 2 W(x,y,z,w), respectively. Before writing down the explicit form of the resulting equations of motion, it is 
useful to summarize the terminology introduced in the course of the above discussion: 

• mean field A (Fig. 3f): denotes the classical soft field as the expectation value of the gauge field A, which is 
induced by the abundance of emitted hard gluons and their collective motion. 

• free propagator iAo (Fig. 3a): refers to the free propagation of hard gluons in the absence of interactions, i.e. 
vanishing coupling g = 0. 

• mean-field propagator iA (Fig. 3b): denotes to the tree-level propagator without quantum corrections, i.e. the 
free propagator with an arbitrary number > of attached external legs coupling to the soft mean field, but 
without closed loops that correspond to quantum self-interactions. 

• full propagator iA (Fig. 3c): terms the dressed propagator of the hard quanta, that is renormalized by both 
the interactions with the soft mean field and the self-interactions among the hard quanta. 

• full vertex functions — igV, —g 2 W (Fig. 3d, 3e): represent the 3-gluon and 4-gluon vertices with the internal 
lines being the full hard propagator including mean-field and quantum interactions. 
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1. Yang-Mills equation for the soft mean field 



The equation of motion for the soft field AJx), is given by @, i.e., STp/SA = -KF> - JC^ o A, ft om which one 
obtains, upon taking into account the initial condition (|2l|), JC^ = 0, the Yang-Mills equation for A: 



(37) 



where = D F — F D with the covariant derivative defined as D X = D X [A] — d x — igA X (x), and F\^ = F\ M L4] = 

[D\ , Dfj\ /{—ig). The second term on the right side is the initial state contribution to the current, according to the 

r A, b. 



condition ©, JC^ o A = Jf p d 4 y /C^f (*, y) A^\y). 
Rewriting the left hand side of ( p7| ) i 



-A, «* -6 



A/ t 



where, upon taking into account the gauge constraint (|^), the —n^n\A in T> Q * A /I does not contribute, because 
= (n ■ A) = n v 'A v , eq. (|37j ) may be expressed in the alternative form (see Fig. 5a): 



{ (V + JC< a >) a}°(x) + S» + ?» = °- 



Here the function 5 contains the soft-field self-coupling, 



S uO) = 3(1) *(x) + 3 (2) „(x) 



(39) 



(40) 



(z) = - 



r 2 — _ 

n d 4 ^ V ^(x, X2) A V ' \x 1 )A X ' c {x 2 ) 



i=l 



%)M= + [ f[d 4 Xi W 0fl t^(x,x u x 2 ,x 3 ) A I/ -\x 1 )A X - c (x 2 )A a ' d (x 3 ) , 



(41) 
(42) 



and the current j is the induced current due to the hard quantum dynamics in the presence of the soft field A: 

= ifl) M (z) + 1(2) ,W + 1(3) ^ ( 43 ) 



1(3) 



- 2 

""I / II^^M^^l,^) A^ fcc (x!,X 2 ) 

" T / II ^ $£(x, x 1 ,X2,x 3 ) A"' \ X1 ) A Act < cd (x 2 ,x 3 ) 

J P i=l 

jf]} ^o^Aa^- zi,*2,z 3 ) A"" > 66 (^,yi) A AA < cc (x 2 ,y 2 ) 



2 JL 

6 ./p 



xA [x 3 ,y 3 ) V 0ll , v , xirr , {yi,y 2 , ys) 



(44) 
(45) 



(46) 



It should be remarked that the function S on the left hand side of (|37]) contains the non-linear self-coupling of the 
soft field A alone, whereas the induced current j on the right hand side is determined by the hard propagator A, 
thereby generating the soft field. 
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a) 



= (D' 1 + % {2) ) A + S + j 



b) 



- + [-1/2 ~< + -« 

X \ X ' vO, X 




+ [-1/2 ^@ + i/2 ig) + 



-i 



o = A' A —(Ag + ^: (2) ) + n + n 



x y 



x y 



x v , < 
+ 17/2 • — • 

5(x,y) 



' x y, 

— i/6 • — •: , 

&(x,y) 



i/2 



— 1/2 



i/6 



$(x,y) 



X s 



— 1/24 




FIG. 5. Diagrammatic representation in terms of the rules of Fig. 3, of the equations of motion: a) The Yang-Mills equation 
(3^), ( |3i|) f or the soft field A with the self-coupling contribution H, eqs. (|io|)-(Jf^) and the generating hard gluon current j, eqs. 
( 200)-( 204). a) The Dyson-Schwinger equation ( [f7|) for the hard propagator A with the mean-field polarization tensor II, eqs. 
( I8|)-(|5o|)~and the quantum contribution II, eqs. ([5l|)-([55|). 



2. Dyson-Schwinger equation for the hard gluon propagator 

From the equation of motion ( p0| ) for the hard propagator, A°5,(i, y), that is, STp/SA = K,^ / '(2i) , one finds after 
incorporating condition j2l|), K,^ — 0, the Dyson-Schwinger equation for A (see Fig. 5b): 

( —l 1 a ^ 

(a) - (Ao)" 1 - JC {2) + n + n (x,y) = 0, (47) 

where A is the fully dressed propagator of the hard quantum fluctuations in the presence of the soft mean field, defined 
by ((2?]), whereas Ao is the free propagator, given by (f24|). The polarization tensor II has been decomposed in two 
parts, a mean-field part II and a quantum fluctuation part II. The mean-field polarization tensor II incorporates the 
local interaction between the hard quanta and the soft mean field, 

EC,0,y) = n (1) a ^(x,?/) + U^^y) (48) 
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n (1) t,(s,y) = f *p(*>v) J p d 4 zV a H x ( X ,y,z)A X ' c (z) (49) 
n (2) t,(^,y) = y J p d 4 zd 4 wW a % d Xr7 (x,y,z,w) A^\z)T'\w) . (50) 

plus terms of order g 3 A which one may safely ignore within the present approximation scheme. The fluctuation 
polarization tensor II contains the quantum self- interaction among the hard quanta in the presence of A. It is given 
by the variation 2iSTp^/SA of the 2-loop part Ip , eq. (|29|), of the effective action Tp, 



ip 



n (2) ^(^,») = - ^f- / II^A, ^^(i.xi.xa) A^'-'^r,^) A-'' dd '(a; 2 ,y 2 ) ^'av(j/2, yi, y) (53) 

i=l 

n (3) t,(^y) = -y j Yldtxidtyi W 0ft f^(x, Xl ,x 2 ,x 3 ) A^'^' ( Xl , yi ) A^' dd ' (x 2 ,y 2 ) 



x A^'-^Cxa.ift,) W^tCift.ift.yi.y) (54) 



n (4) t,( a; 'y) = -^ ^Jldtxidtvidtzi W 0l ^(x >Xl ,x 9 ,x 3 ) A ap '' df ' (x 2) z 2 ) A Tp "' ef " (x 3 ,z 3 ) 



xV p C p f ,' p f (z 3 ,z 2 , Zl ) A px ' >f c ' {zx,yx) A Xa ' :Cd ' (xi,y 2 ) v£*{ yi ,y 2 ,y). (55) 



Note that the usual Dyson-Schwinger equation in vacuum is contained in ( |47| ) ~([55"|) as the special case when the mean 
field vanishes, A(x) — 0, and initial state correlations are absent, K,( 2 \x,y) = 0. In this case, the propagator becomes 
the usual vacuum propagator, since the mean-field contribution II is identically zero, and the quantum part II reduces 
to the vacuum contribution. 



III. TRANSITION TO QUANTUM KINETICS 



The equations of motion (^) or ([39]) for F or A^, and ( fi"7| ) for A^, are non- linear integro-differential equations 
and clearly not solvable in all their generality. However, the field-equations of motion fl37|) or ( f47| ) can be cast into 
much simpler quantum-kinetic equations with the help of the Wigner-function technique and gradient expansion, and 
the assumption of 2-scale separation. As a result one obtains finally the three master equations mcntioncned in Sec. 
1: a simplified Yang-Mills equation decribing the space-time change of A, and two equations for the gluon propagator 
A, namely, first, an evolution equation for the QCD evolution in momentum space, and second, a transport equation 
for the space-time development in the presence of A. In order to achieve this result, one needs to make a third key 
approximation (in addition to the two approximations of Sec. II A), namely, 

• approximation 3: It is assumed that the induced soft field A M is slowly varying on the scale of the short-range, 
hard quantum fluctuations, that is, the gradient of the soft field is small compared to the Compton wavelength 
of the hard quanta. Then one can treat the quantum fluctuations of A(r, k) at short distances separately from 
the collective effects represented by to the soft field A(r) with long wavelength. 
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A. Quantum and kinetic space-time regimes 



The key to derive from (|37j) or ( [47j ) the corresponding approximate quantum-kinetic equations is the separability of 
hard and soft dynamics in terms of the space-time scale A = where /i is the parametric momentum scale introduced 
in (^) . This implies that one may characterize the dynamical evolution of the gluon system by a short-range quantum 
scale r qua <C A, and a comparably long-range kinetic scale r^m A. Low-momentum collective excitations that may 
develop at the particular momentum scale gfj, are thus well separated from the typical hard gluon momenta k± > /i, 
if g <C 1. Therefore, collectivity can arise, because the wavelength of the soft oscillations ~ l/<7£t 1S much larger than 
the typical extention of the hard quantum fluctuations ~ I emphasize that this notion of two characteristic scales 
is not just an academic construction, but rather is a typical property of quantum field theory. A simple example is 
a freely propagating electron: In this case, the quantum scale is given its the Compton wavelength ~ l/m e in the 
rcstframe of the charge, and measures the size of the radiative vacuum polarization cloud around the bare charge. 
The kinetic scale, on the other hand, is determined by the mean-free-path of the charge, which is infinite in vacuum, 
and in medium is inversely proportional to the local charge density times the interaction cross-section, ~ l/(n e o~i n t). 
Adopting this notion to the present case of gluon dynamics, I define r qua and r^i n as follows: 

• The quantum scale r qua measures the spatial extension of quantum fluctuations associated with virtual and 
real radiative emission and re-absorption off a given hard gluon, described by the hard propagator A. It can thus 
be interpreted as its Compton wavelength, corresponding to the typical transverse extension of the fluctuations 
and thus inversely proportional to the average transverse momentum, 

r qU a = A ~ 1 , ( fc_L ) > V , (56) 



where the second relation is imposed by means of the definition (|6|) of hard and soft modes. In general, A 
can be space-time dependent quantity, because the magnitude of (k±) is determined by both the radiative 
self-interactions of the hard gluons and ther interactions with the soft field. 

• The kinetic scale r^in measures the range of the long-wavelength correlations, described by the soft mean-field 
A, and may be parametrized in terms of the average transverse wavelength of soft modes (<7_l), such that 

run = A ~ 7 — r , (q±) < gfi, (57) 

where A may vary from one space-time point to another, because the population of soft modes A(q) is determined 
locally by the hard current j with dominant contribution from gluons with transverse momentum ~ fi. 

The above classification of quantum- (kinetic-) scales specifies in space-time the relevant regime for the hard (soft) 
dynamics, so that the separability of the two scales r qua and ruin imposes the following condition on the relation 
between space-time and momentum: 

A < A, or ( k± > > n > gn w ( q ± ) . (58) 

The physical interpretation of (|5S| ) is simple: At short distances r qua <C l/(gfx) a hard gluon can be considered 
as an incoherent quantum which emits and partly reabsorbs daughter gluons, corresponding to the combination of 
real brems-strahlung and virtual radiative fluctuatiuons. Only a hard probe with a short wavelength A < r qua can 
resolve this quantum dynamics. On the other hand, at larger distances ruin ~ l/(.9/ i ) ; a gluon appears as a coherent 
quasi-particle, that is, as an extended object with a changing transverse size corresponding to the extent of its intrinic 
quantum fluctuations. This dynamical substructure is however not resolvable by long-wavelength modes A > Thin 

of 

the soft field A. 

Accordingly, one may classify the quantum and kinetic regimes, respectively, by associating with two distinct 
space-time points x M and y M the following characteristic scales: 

W = \ (a£ - - gv 

0£ = d£ + W ~ /i . (59) 



(x» +y») ~ A 



9V 
1 
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On the kinetic scale the effect of the soft field modes of A on the hard quanta involves the coupling gA to the hard 
propapgator and is of the order of the soft wavelength A = l/(g/j,), so that one may characterize the soft field strength 
by 

gAp. ( r ) ~ gn, gF^{r) ~ sV , (60) 

plus corrections of order g 2 /! 2 and <7 3 /x 3 , respectively, which are assumed to be small. 
On the quantum scale, on the other hand, 

~ k 2 ± > M 2 » <?V ~ gF^, (6i) 

and one expects that that the short-distance fluctuations corresponding to emission and reabsorption of gluons with 
momenta k± > /i, are little affected by the long-range, soft mean field, because the color force ~ gF acting on a gluon 
with momentum k± ~ ji produces only a very small change in its momentum. 



B. The kinetic approximation 



The realization of the two space-time scales, short-distance quantum and quasi-classical kinetic, allows to reformulate 
the quantum field-theoretical problem as a relativistic many-body problem within kinetic theory. The key element is 
to establish the connection between the preceding description in terms of Green functions and a probabilistic kinetic 
description in terms of of so-called Wigner functions [pl[ . Whereas the 2-point functions, such as the propagator 
or the polarization tensor, depend on two separate space-time points x and y, their Wigner transforms utilizes a 
mixed space-time/momentum representation, which is particularly convenient for implementing the assumption of 
separated quantum and kinetic scales, i.e., that the long- wavelength field A is slowly varying in space-time on the 
scale of short-range quantum fluctuations. Moreover, the trace of the Wigner-transformed propagator is the quantum 
analogue of the single particle phase-space distribution of gluons, and therefore provides the basic quantity to make 
contact with kinetic theory of multi-particle dynamics |18| . 

In terms of the center-of-mass coordinate, r = |(.t + y), and relative coordinate s — x — y, of two space-time points 



x and y, eq. (59), one can express any 2-point function Q(x,y), such as A, II, in terms of these coordinates, 

Qt^v) = e£(r + l,r-l) = g*(r,8), . (62) 

The Wigner transform Q{r, k) is then defined as the Fourier transform with respect to the relative coordinate s, being 
the canonical conjugate to the momentum k. In general, the necessary preservation of local gauge symmetry requires 
a careful definition that obeys the gauge transformation properties 0, but for the specific choice of gauge (||), the 
Wigner transform is simply J32J : 

G(r,s) = J (0e- lfcs G(r,k) , Q(r,k) = J d 4 se lk - s G(r,s) . (63) 

The Wigner representation ( |63|) will facilitate a systematic identification of the dominant contributions of the soft 
field A to the hard propagator A, a concept that has been developed by Blaizot and Iancu 0: First, one expands 
both A and A = A.^ + ^^\A] m terms of gradients of the long-range variation with the kinetic scale r, and second, 
one makes an additional expansion in powers of the soft field A and of the induced perturbation <^A^ ~ 9^-[o]- 



1. Gradient expansion 



To proceed, recall that the coordinate r M describes the kinetic space-time dependence ~ r^m, whereas s measures 
the quantum space-time distance ~ r qua . In translational invariant situations, e.g. in vacuum or thermal equilibrium, 
Q(r,s) in ( p3| ) is independent of r M and sharply peaked about s M = 0. Here the range of the variation is fixed by 
A = 1/fJ,, cq. (|56|), corresponding to the confinement length const, x 1/A in the case of vacuum, or to the thermal 
wavelength const, x 1/T in equilibrium. On the other hand, in the presence of a slowly varying soft field A with a 
wavelength A = l/(g/i), eq. (57), the s^ dependence is little affected, while the acquired r M dependence will have 



a long-wavelength variation. In view of the estimates (|59|), one may therefore neglect the derivatives of G(r,k) with 
respect to r M which are of order <?/i, relative to those with respect to s^ which are of order /i. 
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Hence one can perform the so-called gradient expansion of the soft field and the hard propagator and polarization 
tensor in terms of gradients (s ■ d r ) n , and keep only terms up to first order n = 1, i.e., 

Mx) = I M (r + |) ~ A,(r) + I ■ d r A»{r) , (64) 

and similarly for A^y) = A^(r — s/2), as well as 

A^(x,y) = A flv (r,s) ~ A^(0,s) + s ■ d r A M „ (r, s) (65) 

Uftv (x, y) = (r, s) ~ (0, s) + s ■ d r (r, s) , (66) 
Then, by using the following conversion rules [0,E!3| to carry out the Wigner transformations, 



d 4 x'f(x,x')g(x',y) => f(r,k)g(r,k) + - {(d k f) ■ (d r g) - (d r f) ■ (d k g)} (67) 

h(x)g(x,y) => h(r)g(r,k) - % - (d r h) ■ (d k g) h(y)g(x,y) h(r)g(r,k) + % - (d r h) ■ (d k g) (68) 

d£f{x,y) => (-ifc" + \d?)f{r,k) d£f{x,y) (+ik^ + \d?)f{r,k) , (69) 
the transformed polarization tensor II (r, k) is obtained from \\{x 7 y), eqs. ( [l8| ) and (|5l|), with 

IV (r, k) = IV (r, k) + IV 0, k) (70) 



where the soft mean-field contribution [c.f. (|49|), (50)] is 

<,(r,fc) = (n (1) + n (2) ) ab (r,k) (71) 

n ( i) ^(r,*) = |<tA(fc,0,-fc)A A ' c (r) (72) 

n (2) °^(r,fc) = - l {w o a % d Xa (k,0A-k)A X ' c (r)T- d (r) , (73) 

and the quantum contribution [c.f. (|5l|)-(|5^)1 is 

n^(r,fc) = (n (1) + n (2) + n (3) + n (4) )° 6 (r,k) (74) 

^ { i)U r ^) = + l -Y J W 0lx l b x i(k,q,-q,-k)A^> cd (r,k) (75) 

n (2) <tM) = + y / ^0^(^-9,-3') A AA '' cc ' (r, g) A aa '' dd ' (r,q') V d !y b (r;q',q,—k) (76) 

■ tt> ;., _ .'/ / " '/ a P w acde/u „ „' „\ A AA'.cc'^, \tr<r',ddf, 

x A^'^ e '(r,p) ^; c A ; fc t/ (^9,'z',P,-fc) (77) 
-<?,-</, -p') A°o'- d f'(r,q) A T P"^"{r,q') 

xV p C p f ' f (r ]q ,q',-p) A> A ^ c >,p)A ACT '' cd '(r,p') (r; p, p', -fc) . (78) 



V>' fe ) = -fi / (2^ < (2$ i W 'o^(*,-ff.-9'.-l/) A ap ' ' df ' (r,q) A Tp " ' ef " (r,q') 



Here the 3- and 4-gluon vertex functions from (g4J), (35) and ( |36|) depend explicitly on r, 
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V{r;h) =V (k i )+O(g i f(r,k i )) W(r-k l ) = W Q (k i ) + 0(g i f(r,k l )) 

with the bare point-like vertices Vq, Wq being r-independent and given by 

V %Aki,k 2 ,k 3 ) = -tf abc { g^lhx - k 2 ) v + g^(k 2 - k 3 ) x + 9l/X (k 3 - k 2 )„) 

W^AklMMM = - { {r e f bde -f ade f Cbe ) 9X»gv a + (f<*° fed* _ fade f bce) 



(79) 
(80) 
(81) 



With the above formulae, one can now convert both the Yang-Mills equation ( |37| ) and the Dyson-Schwinger equation 
( p7| ) into a set of much simpler equations. For the Dyson-Schwinger equation, the Wigner transformation together 
with the gradient expansion yields two distinct equations for the hard propagator A A " y (r, k), namely, (i) an evolution 
equation [j, and (ii) a transport equation. They are obtained fll],H by taking the sum and difference of Wigner- 
transform of ( |47| ) and its adjoint, using the rules (|67|)- 
(i) evolution equation: 



(ii) [k'--dt] A^(r,k) 



- {n£, A™}(r,k) + 



(r, k) 



i 
4 



(r, k) 



a A i> , d k x A° 



(r,k) (82) 



(ii) transport equation: 



(k-d r ) A^(r,k) + - 



K , A™ 



n£ , a° v 



(r,k) - {d*K,d h x A™}(r,k) 
(r,k) + \ {d£nZ,d r x A«»}(r,k) - \ {a r A n^,^A^}(r,fc) 



(83) 



where d 2 r = d r ■ d r , [A,B] = AB - BA, {A,B} = AB + BA. In (||) and (|3|), l P = 1 (0) for A F , A F (A>, A<) 
arises as the transform of Sp(x, y). The function d^ u (k) is the sum over the gluon polarizations s, (emerging from the 
Fourier transform of the operator (p5|)), 

(84) 



d^k) = £ e,{k,s).el{k,s) = n ^+^ + {n 2 +ak 2 } Jh^ 

t—' n ■ k (n ■ ky 



8=1,2 



with the properties d^(k) = 2 Fl k^d^Qt) — ► and n^d^ v = = d^ v n v . Furthermore, the initial state contribution 
appearing in (|37|) and (|47|), which contributes only at r° — to, has been absorbed into the hard propagator, 



A-}(r,k) ee A-J(r,k) - £$(r, k ) 5(r° - to) 



(85) 



For the Yang-Mills equation (pTj), determining F (r), one obtains on the same level of approximation a compact 
expression in terms of the hard current j Pc 



D X ' ab ,F" Xtl 



(r) = -J» 
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D x ,A„ a (r, k) 



(86) 



where 7 m „a<x = tg^gxa - g^xg V a - g^ <T g v \®n&D x = d r x -igA x {f). 



t In Ref. |l^] the 'evolution' equation was termed 'renormalization' equation, a term that may be misleading. In order to 
avoid confusion with the 'renormalization group' equation, the name evolution equation appears more suitable. 

* This property reflects that in the non-covariant gauges (C6) only the two physical polarization states propagate, i.e. those 
with with Epfc^ = 0. For comparison, in the covariant Feynman gauge, d M " = g^ v , = 4, and k^d^ = k v 7^ 0. 

§ Note that in the kinetic approximation, the piece eq. (kq), does not contribute, because it has two additional A 
insertions and is down by a factor g/^ 4 as compared to and 
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2. Expansion in powers of gA 



In order to isolate the leading effects of the soft mean field A on the hard quantum propagator A, I follow Ref. 
to separate the quantum contribution from the mean field contribution on the basis of the assumption that the A 
field is slowly varying on the short-range scale of the quantum fluctuations. To do so, recall eq. (|27|), 



A(r,k) = A@(r,fc) + 6A m (r,k) 
with the quantum piece A^ and the mean-field part defined by 



V = A " 



A=0 



= A - 



n 



A=0 



6A m = A 



n 



(87) 



(88) 



and the free-field propagator Ao and the mean- field proagator A are given by (J2J) and (|26|), respectively. Given the 
ansatz (jS7f), with the feedback of the induced soft field to the hard propagator being contained in 5 A ^ , the latter is 

now be expanded in powers of the soft field coupling gA, and it is anticipated that the mean-field induced part 8A^, 

is a correction being at most g times the quantum piece A^, that is, 



<5A [3] (r,fc) = ]T ± (gA(r) ■ d k ) n A [Tj] (k) ~ gA(r) - 9 fe A p] (r, k) , 



and, to the same order of approximation, 

d^A^^k) ~ ff (^)^A m/w (r,fc), 



(89) 



(90) 



where, on the right side, the space-time derivative acts only on A. Now the decomposition (|87j ) with the approximation 
( [sit ) is inserted into eqs. (|82]), (|8^), (ftq) , and all terms up to order g 2 ^i 2 A^ are kept. The resulting equations can 
be compactly expressed in terms of the kinetic momentum rather than the canonical momentum k^ (as always 
in the context of interactions with a gauge field psfi ), which for the class of axial- type gauges (||) amounts to the 
replacements 



kn 



K„ 



9 Mr) , 



D 



dl gd r JF(r)d k v 



—2 ■ 



(91) 



Taking into account the approximation 3 of Sec. Ill A implying K 2 A 3> D" r A, one finds for the the evolution- 
transport-, and Yang-Mills equation, eqs. (p2h, (§3H, and (p6|), respectively, 



{ X 2 , Ag} (r, X) = ^(#) + i {n^, A^} (r, K) 

a 



if • D r , A^ 



(r,*T) 



(92) 

2 L" CT ' J [o]J v, " ; 2 L" CT ' "Pi J v ''" ; ( ' 93 ' ) 
to =-%(r) = -g j 0^ Tr{ -^A„ "(r.iiQ + A^ "(r,tf)^} , (94) 
where the color indices are suppressed, noting that A a ^ v = A M „, Fx^ = T a F^ M , j M = T a j®, and Tr[. . .] = TrT a [. . .]. 



D r , F\fj. 



3. The physical representation 

One sees that the original Dyson-Schwinger equation (f47p reduces in the kinetic approximation to the set of algebraic 
equations ( |92"|) and (|93|). Now recall (c.f. Appendix B3)that in the CTP framework these equations are still 2x2 
matrix equations which mix the four different components of A = (A F , A>, A<, A F ) and of II = (U F ,U> ,U< ,Il F ). 
For the following it is more convenient to employ instead an equivalent set of independent functions, namely, the 
retarded and advanced functions A ret , A adv , plus the correlation function A cor , and analogously for II. This latter set 
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is more directly connected with physical, observable quantities, and is commonly referred to as physical representation 
R: 



r adv 



(95) 



Similarly, for the polarization tensor the retarded, advanced and correlation functions are defined as (note the subtle 
difference to (95h): 



jjrei = U l 



IT 



if 



jjcor = _ n > jj 



(96) 



Loosely speaking, the retarded and advanced functions characterize the intrinsic quantum nature of a 'dressed' gluon, 
describing its substructural state of emitted and reabsorbed gluons, whereas the correlation function describes the 
kinetic correlations among different such 'dressed' gluons. The great advantage jl^j23] of this physical representation 
is that in general the dependence on the phase-space occupation of gluon states (the local density) is essentially 
carried by the correlation functions A > , A < , whereas the dependence of the retarded and advanced functions, A ret , 
A adv , on the local density is weak. More precisely, the retarded and advanced propagators and the imaginary parts 
of the self-energies embody the renormalization effects and dissipative quantum dynamics that is associated with 
short-distance emission and absorption of quantum fluctuations, whereas the correlation function contains both the 
effect of interactions with the soft mean field and of statistical binary scatterings among the hard gluons. 

In going over to the physical representation, one finds then that eqs. (^||) and (|9^) give a set of 'self-contained' 
equations for the retarded and advanced functions alone, 



L J [iv 



1 / ret **> ret ^ ret ret 

— ^ IJgm2u ^adv -|- ^adv Yl a dv 



K ■ D r , Aadv 



= — — ( fjadv A adv 
[iv 2 



plus a set of 'mixed' equations for the correlation functions 

1 
2 



n< a q 



dv 



rr e * a^ 



K D r , A< 



= - - ( IT 4 A adv + n ret A^ 



^ret ret 
/\ adv Yladv 



jjadv 
jradv 



fjLV 

A ret 

- A ret U^ 



(97) 
(98) 

(99) 
(100) 



The equations (|97|)-( 100 ) may be further manipulated by the following trick. Let the imaginary and real components 
of the retarded and advanced propagators be denoted by 

1 



2ImA„ 



A a 



dv 



[IV 



ReA^ 



A a ' 



dv 



fl.V 



(101) 



with A ret = (A adv )* and 6(K°) A ret — 9(—k°)A adv . The analogous decomposition of the polarization tensor in 
terms of its real and imaginary components defines the quantum part II as the sum and difference of the retarded 
and advanced contributions, respectively, 



Tau = 2ImLV = i n ret -n 



)J,V 



RelV = I (n ret + fl adv ) 

2 V / llv 



(102) 



and similarly for the mean-field part II, associated with presence of soft field. The imaginary parts p and T are the 
spectral density and spectral width, respectively, of the hard gluons. 

In terms of this representation one obtains from eqs. (|93|)-(|9^) and (97)- (100) the following final set of master 
equations: 



{ K\ p}^ = {Ren, p] ^ + {f , RcAp,} ^ + g (f^ p Xv + p* F Xl 



= +i 



IT 



ReA 



[o] 



9 K X F 'd*^ 



z[Ren,A- 
■ 9{fIa^ Xv 



^{ ficor 'C-^{ f '^rL' 



A c °? x F Xl 



[o] , 



D r , F Xfl 



= -3k = 



Tr 



{ (-K»A< 



A c ° r v K, 

H 1 



-)} 



(103) 

(104) 
(105) 
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The physical significance of the ( |103| ) and (104) is the following Eq. ( Jl 03] ) determines, in terms of the spectral 
density p, the state of a single gluon with respect to its virtual fluctuations and real emission (absorption) processes, 
corresponding to the real and imaginary parts of the retarded and advanced polarization tensor in the presence of the 
soft field F. Eq. (104), on the other hand characterizes, in terms of the correlation function A cor , the correlations 
among different such gluon states. The polarization tensor appears here in distinct ways. The first two terms on 
the right hand side account for scatterings between the single-gluon states. The next two terms incorporate the 
renormalization effects which result from the fact that the gluons between collisions do not behave as free particles, 
but change their dynamical structure due to virtual fluctuations, as well as real emis sion and absorption of quanta. 
The last two terms account for the soft interaction with the mean field F. Eq. (105) finally determines the rate of 
change of he soft field F by the hard gluon current, which involves the full correlation function A cor . 

The interlinked structure of eqs. (103 )-(|l0"lj) is very convenient for explicit calculations (demonstrated in Sec. 4). 



It provides a systematic solution scheme, as discussed below, to solve for the three quantities of interest, namely, the 
spectral density p, the correlation function A c 



and the mean field A. In view of (h03j)-(105) the natural logic is a 



stepwise determination of p 



A cor 

A [o] 



5A™ r 

[A] 



J 



C. General solution scheme 



Let me exemplify the above interpretation of fllOSj ) and (104) in more quantitative detail (see also Refs. |19| , p2| ). 
The formal solution of ( 103 ) for the retarded and advanced functions is |19| , 



a ret A re ^ 

[iv p.f 



(a™ 1 it 



/(/(• 



\adv A aQ 

fiv fiv 



dv a &dv 



(106) 



_ret 
jadv 



where Uadv = Uadv -f n . This determines p^ LU via (101). Once Aadv is known, the solution of (104) for the 
correlation function is given by |l9| , 



a cor 

^■fJ,V 



ret /\ cor — 1 ^adv 



A ret A 



A re * JJ cor /^adv 



(107) 



with n cor = U cor + n .It has the general form [g2j 

A c ^(r,K) = -ip^(r,K)G{r,K) , 
i.e., the convolution of the spectral density p^ v with the phase-space density of hard gluons G, 

G(r,K) = 1 + 2 S (r,K) 



(108) 



(109) 



where the 1 comes from the vacuum contribution of a single gluon state, and the 2g represents the correlations with 
other hard gluons that are close by in phas e-space. Note that the function g is constrained to be a real and even 
function in K (c.f. Appendix F). From (108) it follows that the total number of gluons N in a space-time element d A r 
is 



d 4 K 



(27T) 



Tr 



d 4 K 

(27T)* 



p(r,K)G(r,K) , 



(110) 



where dfj, v (K) is the polarization sum given by (p4]), d^ 1 = 2d fJ-u , and p — ^d^p^ and an averaging over the 
transverse polarizations and the color degrees of freedom is understood. 

The above formulae become immediately familiar, when considering for illustration the simplest case of a non- 
interacting system of gluons, the free-field case. In this case, II = and one finds, utilizing the formulae of Appendix 
F, for the free retarded and advanced functions: 



A ret I ts\ (K) 



A adv I d^uiK) 



K 2 -if.' 

Hence, the free-field spectral density po which is the difference between Aq 6 * and AQ dv , is on-shell, 



(111) 
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-ip 0fiv = A r % v - Ag* = 2ir6(K 2 )d^(K) , (112) 

by means of the principal- value form ula (K 2 ± ie)^ 1 = PV(l/K 2 ) =F iirS(K 2 ). The free- field correlation function Aq 01 " 
is then readily determined via ( |l08|) , 

A c °^(r,K) = -2m6(K 2 ) G (r,K) d M „(K) , (113) 
and so, with Go = 1 + 2go, the number of on-shell gluons per d A r is 

Mr) = J { J" 3 K 2K o G o(r,K) , G (r,K) = G (r, K)\ K0=]jt] = (2^) 3 2X°|| (114) 



The free-field exercise, eq 



acting case, eqs. (|109)-(110) 



s. (111)-(114), illustrates the two main properties, which hold also for the general inter- 



(i) The spectral density p^(r, K) describes the 'dressing' of a single gluon state with momentum K with respect 

to its radiative quantum fluctuations, i.e., its fluctuating coat of emitted and reabsorbed gluons. The function 
Tr[d~^p fiv ] is the intrinsic gluon distribution, that is, the number of gluons inside this gluon state. The spectral 
density is a property of the state itself and therefore is nonvanishing even in vacuum, in the absence of a medium. 
For on-shell particles, p^ v oc 5{K ) and therefore there are no intrinsic gluons present. 

(ii) The correlation function A^fJ(r, K) describes an interacting ensemble of such fluctuating gluon states, and is 

given by the number density G(r, K) of those gluons weighted with their spectral density /9 M „, containing the 
intrinsic gluon density of each of them. For the non-interacting case, it obviously reduces to an ensemble of 
on-shell particles with K° = \K\. 

In closure of this Section, a generic solution scheme may be the following iteration recipe (which is exemplified in 
the next Section): 

1. Solve the evolution equation ( p_03[ ) for An* and the associated spectral density p at starting point t = to with 
specified initial condition p(to) — po at a large initial momentum or energy scale Q. This can be done just as in 
free space, except that the kinetic momentum K = k — gA carries now an implicit dependence on the soft field 
A with specified initial value A(to). 



2. Solve the transport equation (104) for the correlation function A cor = A^ + SAg?. This involves, a) the 



construction of A^ r with the help of p and Aadv from step 1, and b) the calculation of the mean-field induced 



correction ^Aj^ r from the right side of (104). The resulting space-time evolution of A cor describes then the 



evolution of the gluon density G within a time interval between to and t\ ~ l/(Kj_i), corresponding to the 
evolution from Q down to a mean Kj_i at t\. 



Insert the solution for the full correlation function A cor into the current j p on the left side of (105) and integrate 
over all momenta K from the initial momentum scale Q down to the hard-soft scale p. This gives the current 
induced by the motion of the total aggregat of hard gluons during the evolution between to and t\. Then solve 



the Yang- Mills equation (105) to determine the soft field A (equivalently F) that is generated at t\ as a result 



of the hard gluon evolution. 
4. Return to 1., and proceed with second iteration, replacing A(to) by A(ti), and so forth. 



IV. SAMPLE CALCULATION: HARD GLUON EVOLUTION WITH SELF-GENERATED SOFT FIELD 



This Section is devoted to exemplify the practical applicability of the developed formalism by following the solution 
scheme of Sec. Ill C for the specific physics scenario advocated in the introduction and schematically illustrated in 
Fig. 1 . I consider a high-energy beam current of hard gluons as it evolves in space-time and momentum space, and 
eventually induces its soft mean field. 
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A. The physics scenario 



(i) The initial state is modeled as an ensemble of a number Ao of uncorrelated hard gluons. The Lorentz-frame of 
reference is the one where the gluons move with the speed of light in +z-direction. The initial gluon beam is 
prepared at 

Tq = (ta,fj_ , z ) , t = z = , < rj_ < R , (115) 



corresponding at to to a sheet located with longitudinal position zo with transverse extent up to a maximum R, 
specified later. 

(ii) The initial hard gluons are imagined to be produced at some very large momentum scale Q 2 3> A 2 , with 
their energies and longitudinal momentum along the z-axis being ~ Q. These gluons are therefore strongly 
concentrated around the lightcone with momenta 

k° ~ k 3 « Q , ft « , (116) 



and hence have very small spatial extent Ar ~ 1/Q. That is, the initial state gluons are taken as bare quanta 
without any radiation field around them. 

(iii) The subsequent time-like evolution of these bare gluons proceeds then by two competing processes: a) the 
regeneration of the radiation field by emission and re-absorpton of virtual quanta, and b) the brems-strahlung 
emission of real gluonic off-spring. As a consequence, phase-space will be populated with progressing time by 
more and more gluons. The typical energies decrease, whereas the average transverse momentum increases (c.f. 
Fig. la), but yet within the hard momentum range 

Q 2 > k\ > (i 2 > A 2 . (117) 



Eventually, the evolving gluon system reaches the point at which the transverse momenta become of the order 
of the energies. This point is defined to be characterized by the scale fi - the transition from hard, perturbative 
to soft, non-perturbative regimes. When K\ < [x 2 , the individual gluons cannot be resolved anymore, and their 
coherent color current acts as the source of the soft mean field. 

(iv) Because of the restricted kinematic region ( |ll7[ ) of the hard gluon dynamics, the coupling a s = <? 2 /47r satisfies 

a s (k 2 ± ) < 1 a s (k 2 ± )ln(Q 2 /k 2 _) ~ 1 for all k\ > fi 2 , (118) 



so that a perturbative evaluation of the hard gluon interactions is applicable, provided [i > 1 GeV. The per- 
turbative analysis in the following subsections will be restricted to leading order: the hard gluon interactions 
then includes only radiative self-interactions ~ g 2 , but no gluon-gluon scatterings ~ g 4 R, or other higher- loop 
contributions. Hence, for the hard gluon propagator A = A^ + <5A^ of (|7|), (|S8"|), the required accuracy for 

the quantum contribution A^ is, 

Ar* = A + a( 5 2 ;A ) + 0{g A ) , (119) 



where Ao is the free-field solution. On top of this the interaction of the hard gluons with the soft field is treated 

~~ ]), to leading order ~ gA to the solution A^ of (119), 



as a correction £>Apj as in 



8 Aft = C(gA;A m ) + 0(g 2 A). 



Ph 



(120) 



** Aside from g 4 <C g 2 , the neglect of scatterings is reasonable here, because for a beam of almost colinearly moving gluons, 
the scattering cross-section o gig2 oc a s \v\ — 1)2 1 M gi , g2 (ki, k%) is negligible for vanishing flux \vi — Wa| ~ 0. 
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Although this so defined physics scenario, with an initial state of bare gluons, being only stastistically correlated 
and incoherent, may appear to be rather academic, it has in fact valuable physical relevance: For example, it may be 
viewed as the idealized version of the initial density of materialized gluons in the very early stage of a high-energy 
collision of two heavy nuclei. In this example, one expects the materialization of a large number Ao of virtual gluons in 
the wavefunctions of the colliding nuclei, to occur very shortly after the nuclear overlap by means of hard scatterings. 
If one imagines the time of nuclear overlap equal to to — 0, and assume the average momentum transfer of initial hard 
scatterings w Q 2 , then the above idealistic scenario acquires a more realistic meaning. 



B. Choice of lightcone gauge and kinematics 



For purpose of calculational convenie nce , I wil l henceforth work in the lightcone gauge which is a special case of the 
axial-type gauges (||). It is defined by (C3)-(^6|) of Appendix C, that is, 



n-A a = 0, n 2 = (Al=Al,a*), 
corresponding to the gauge fixing term in (|^), 

1 



(121) 



I GF [n-A] = J p d 4 x (~[r 



■A a {x)\ 



with i 



. 



IP \ 2a 

I choose the lightlike vector n M parallel to the direction of motion of the gluon beam along the forward lightcone 

r/ 1 = (n°,n_L,n 3 ) = (l,0j_,-l) 
and employ lightcone variables, i.e., for any four- vector 

— (v + , v~ , v±) 



(122) 



2 + - 

V — V V 



v± 



(v\v 2 ) 



v± 



1 



v^w^ = — \v ' w + v W ) — Vj_ ■ Wj_ 



Then = (n + ,n ,n±) = (0,1, 0j_), so that the gauge constraint (121) reads 

n ■ A = A + = A- = , 



and the non- vanishing components of the gauge-field tensor F 



p u 



-b are 



F + = -d+A 



F +l = d + A l 



= d~A -d l A -ig 



A ,A 



A\A ] 



F = d l A - d°A - ig 
where 9 ± = d/d^ and the index i = 1,2 labels the transverse components. 



Finally, the kinematic imposition (116) reads in terms of lightcone variables 

K+R- = K 2 + K 2 , < (K+) 2 



(123) 



(124) 
(125) 
(126) 

(127) 



(128) 



(129) 



K+ ~ 2AT° ~ 2 AT 3 K~ 







K\ » K 2 



(130) 



Physically this implies that the hard gluons are effectively on mass shell, i.e., their actual virtuality (degree of off- 
shellness) K 2 is small compared to Kj_, the transverse momentum squared, and negligibly small compared to the 
scale (K + ) 2 . Within this kinematic regime, I henceforth consider K 2 /(K + ) 2 — > 0. 
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C. Properties of A M „ and n M „ in lightcone representation 

The most general Lorentz decomposition of the polarization tensor II = II + II in lightcone gauge can be written 
as Uf u (r,K) = S ab U^(r,K), with 

„ , ^ ( K U K U \ „ (K U K V \ „ f n u K v + K u n v \ „ / K 2 n u n u \ „ 

^(r,K) = [g^—lpr) n, + (-^J n u + ( ^ j Hl + ((^#J 112 (m) 

where II||, III, II2 are scalar functions of dimension mass squared and depend on the four- vectors r M and — 
— gA^ . In lightcone gauge, the Ward identity for the gluon propagator |27]] 

lira A-(n-K)n»A^ + K v \ = (132) 

enforces II M „ to be transverse with respect to and symmetric in its arguments and indices, 

n"Il£ = = n>" njt = Kl , (133) 

which implies that 

n,| = -iii = +n 2 . (134) 

Therefore, with n ■ K = K + , 

*br V \ _ X ab ( K/mKu \ tt , cab I K n K v , n n K i> + K^Uu , ^ 



n- (r, iO = ^ - ^ ) Ex + + + ^= ) n„ , I :$r, , 



1 / 



2 



n± = g U« + 7^+)2 n » n » ) nMV n n = sn ± - . (136) 

The corresponding full gluon propagator is given by the inverse of (A) -1 = (Ao) -1 — II. Using the free-field form 

A^(r,K) = S« b d,AK) A (r,K) , d, u (K) = g, v n » K »+ + K » n " + , (137) 



with the scalar functions A = (Ag e *, Ag d ", A c or ) [c.f. eqs.(lll)-(fil|)] 



= K^±ie ' A §°'>,^) = -2mS(K 2 )G (r 7 K) , (138) 



finds 



A ^(r,K)-5 A (r,K) ^_ n±/ ^j ^ + U - (n± - Hn)/^ J f " ( 139 ) 



Now, because of (13C), the last term in ( |139| ) vanishes for K /{K + ) — > 0, and the full propagator A can be expressed 

as the free-field counterparts Ao times a scalar formfactor function Z whose momentum dependence contains only 
the Lorentz invariants n ■ K = K + and K\ : 

Af v {r,K) = 6 ab d^(K)Ao(r,K) Z(r,K+,Kl) , (140) 

where, because of K 2 /(K + ) 2 — > 0, the function d^ v reduces now to 

, , n,,K v + K n n v 

= 9fiv - k+ (141) 

and the formfactor Z is related to the polarization tensor by 



2G 



Z(r,K+,Kl) = 



1 

1 - U±/K 2 



with boundary condition 



Z{Q,K + ,K\)\ 



K=Q 



1 



(142) 



(143) 



Here Q is the renormalization point, determined by the momentum scale of the initial state hard gluons (which 
specified in the next subsection). 

The great advantage o f the l ightc onc gauge becomes evident now: the solution of the full retarded, advanced 
and correlation functions ( 106 )-( 108) boils down to calculating a si n gle sc alar function for each of them, namely the 
formactor Z, which is simply multiplied to the free-field forms (111)-(113). For the retarded and advanced functions, 
with 



one has 



ret 
£ adv 



1 



1 - W^ v /K 2 



K 2 + te v ' ' ± ' u ^ v ' K 2 -ie 

which satisfy the useful relations A r $ = (A£f )* and A r ^(K°,K) = Afj{-K , K). Defining 

— ^ret ^adv 

the spectral density follows immediately as 



(r,K) = d^(K) (-2iri) ZP(r,K+,K 2 ± ) 



p^(r,K) = i(A ret -A adv 

and the correlation function is obtained as 

A™{r,K) = d^K) (-27rt) Z"(r,K+,Kl) (l + 2g(r,K)) . 



(144) 



(145) 



(146) 



(147) 



(148) 



D. Specifying the initial state 



To fix the initial conditions for the scenario described in Sec. VI A, both A and A have to be provided with initial 
values at r° = to = 0. The initial condition for the hard propagator is chosen as 



A „ v (0,r x ,K) , 



(149) 



referring to a statistical ensemble of bare gluon state s at t ime r° — 0, which can be characterized by a single-particle 
density matrix of the Gaussian form as given by eq. (B14) of Appendix B. This ansatz corresponds to an initial state 
source term in (fi"7j) of the form 



K)\ 



(r=(0,?_i_,0) 



JC$(r,K)6(r°)6(r 3 ) = Po ^{K) G (r, K) 



(150) 



As assumed in Sec. IV A, the initial ensemble consists of a total number Mo of bare gluons with total invariant mass 
Q 2 , all moving with equal fractions of the total momentum Q} L = /Mo- That is, each gluon moves initially with 
momentum — (Q + ,0,0jJ colinearly to the others along the lightcone. Throughout the ultra-relativistic limit is 



understood, "Q J 



i.e. Q + ^ A, where A 0.2 — 0.3 GeV. The spatial distribution of these Mq initial gluons 



at r u = is taken as a (^-distribution along the lightcone at 



,3 = 



zq = 0, and a random distribution transverse to 



the lightcone motion. That is, the initial multi-gluon ensemble is prepared at lightcone position r~* 



ir 



0, 
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lightcone time r = to — z o with a transverse smearing < rj_ < Nq / \[Q 2 ^ where th e ty pical transverse extent of 
each gluon is Sr± ps 1/Q 2 *C 1 fm. Accordingly, the initial state spectral density pa in (15C) is taken as 



(if,) 



dK+dK-d 2 K ± 



po(if) = 1 ■ (151) 



The corresponding retarded and advanced functions Aq are of the form (111) 



A adv {K) = pv| _ J zp _ po(X) 



(152) 



Finally, the initial state correlation function Ag or is the convolution of po with the density of bare gluons at the scale 
Q and lightcone time (position) r~ = r + = 0, 



A^,(r,X) = / d 4 A"(2^) 4 V(K') p (K') G (r,K) , 



(153) 



where G (r, K) = G (r)G (K) with 



GoW = — *(r-) <5(r+) (l -^ig 2 ) 



G (X) = = 



(27T) 



(5 - Q+) <5 if - - 



K+ 



S 2 K X 



(154) 



The visualization of the initial gluon density Go in (154) is a 2- dimensional color-charge density: It is spread out in 
the two transversee directions f± in a disc with radius R = 1/s/AfoQ 2 — 1/Qi an d a delta- function in longitudinal 
direction at r + = at time r~ = 0. The normalization is such that the total number Ao of initial bare gluons is given 
by 



dr-dr+d 2 r_ L G {r,K) = Af Q G Q {K) 



(155) 



Finally, because of this statistical ensemble of almost pointlike, bare gluons, one does not expect any collective 
mean- field behaviour at initial time r° = to = and at large Q 2 , so that the magnitude of the soft field is initially 
equal to zero which is consistent with (|12|): 



,r±)\ r —_ 



-0=0 



(156) 



This co mplet e s the construction of the initial state, starting from which, I now address the solution of the set of 
equations 



E. Solving for the spectral density p M „ 



To find the spectral density p^ v , the solution of A^j* and Ajgj" is needed. 



The first correction to the free-field 

solution (112) arises from two contributions: a) from one- loop hard gluon self-interaction of order g 2 that is contained 



in the hard polarization tensor II, and b) from the coupling of the hard gluon propagator to the soft field A in II 
which is of order g A. Within the perturbative scheme (119) and (120) the retarded and advanced propagators are to 
be evaluated to order g 2 from eq. (106) with the internal propagators in Hadv taken as the free-field solutions, 



^ ret 
A a>dv 

A [0] 



ret 

Afl dv 



+ A$ dv ffadv[g 2 ; A ] Ag dv , 



(157) 



^ ret 

with the subsidary condition [K ■ D r ,A^ v ] = + 0(g A ). To order g 2 , the gluon polarization tensor II as given 
by (|74|)-(|7§|), reduces to the one-loop term IT 2 ), because the tadpole term IT 1 ) vanishes as usual in the context of 



28 



dimensional regularization f27|] ; and the two-loop terms H( 3 \ n' 4 ' are of order <? 4 . Hence, (ll ret — H adv ^J ^g 2 - Aq] ii 
(157) reduces to 



(fl ret -U. adv Y (r,K) 

V / LIV 



d 4 q 



V « X C *(K, -q, -K + q) V ( t c > b yjr; K-q,q, -K) 



(2tt) 4 ' UaiAo 

5 cc 'S dd ' d xx '(q) d™'(K-q) { A^(r,q) A c ° r (r,K-q) - A%> r (r,q) A r et (r, K - q)} , (158) 



where the Af (r, if) = A (r, K) d^ v (K) 5 ab are the zeroth order solutions fllll|) and (113) 



The mean-field contribution (|7lj)-(|73j) to the retarded and advanced components of II, on the other hand, vanishes, 
because F ^ = T°F^„ is antismmctric and traceless 0. 



ret ^ ret \ au ret /I _ 

n'"'''A"-/' ] (r,K) = -2gS ab d^ A§ dv (r,K) i^-g px F pX (r) 



. 



Hence, the dependence on the soft field F \ v or is resident only implicitely in the kinetic momentum 



(159) 



so that eq. ( |103] ) becomes formally identical to the case of A — 0, in which — Exploiting this formal analogy, 
one can evaluate explicitly U r u e J — IK^" in the kinematic range Q 2 > (K + ) 2 ^S> K\ > /i 2 by using standard techniques 



of QCD evolution calculus |ll|35[|. Inserting into (158) the free-field expressions for AQ et , Aq , and Ag or , from eqs. 

ret 



(111), (113), one finds that to 0(g 2 ) the polarization tensor Uadv does not depend on r, hence one may write 



IladvfaK) =ILadv(K) Zadv( r ,K) = Zadv(K) 



Using the lightcone variables (126), for the momenta, together with the lightcone phase-space element 



4 "g 2 ^0" 1 + _ 2 sf + - 2 \ n d 1 + j 2 



2 ^ 



(160) 



(161) 



and using (p6|) and (|l42|) , n™* -Hf 11 = Kllodo)^ one finds the formfactors Zcwfa = (1 - n^/if 2 )- 1 ) to leading-log 
accuracy: 



Z ret (K+,K 2 j_) 



1 



exp ■ 



K +2 =Q 2 



dq\ 



K 2 



dq+ a s (q 2 ± ) 



q+ 2nqj_ ' \K+ 



(162) 



where 



Z adv (K+,Kl) 9(K^ 



Z ret {-K+,Ki) 6{-K+) , 



a s (ql) 



4tt 



7 (z) = 2N C z(l-z) + 



1 - z 



(163) 



(164) 



llln^/A 2 ) " V'" v " '"' ' * '1-2, 

and z = q + / ' K + , 1 — z = q + / K + , g = K — q. The effective formfactor function iJ p can be approximately evaluated, 



Z p (K + ,K 2 j_) = Z r 



Z a 



dv 



for Kl> nQ 



Iln 2 (f)-ln 2 (fl)]} for K\ < nQ 
Substituting Z p into eq. ( |l45| ) for A™* and A^ 11 , one obtains for the spectral density p — id~£ ^A^ 



(165) 



et A adv 



^ Note, however, that this cancellation occurs only in the lightcone gau ge (|l22[) with gauge parameter a = 0. In a general 
non-covariant gauge with one encounters on the right hand side of (|159|) a finite term a (n ■ d r ) (n^A^r) + n M j4.„(r)). 
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p[K + ,Kl)=Z"{K+,Kl) 



K+ 



6 (K+ - Q+) S (Kj 



dq 2 ± a s (ql) 



A' 2 



<4 



2tt 



<iz 7(2) p 



A+ 

2: 



(166) 



The previously advocated interpretation of the spectral density jo of an initial state gluon as the density of its 'intrinsic' 
gluon fluctuations becomes clearer now: p represents the structure function of a gluon that was initialized as a bare 
state at Q 2 . Looking at this gluon sta te w ith a resolution scale Kj_, o ne se es at K\ = Q 2 only the initial bare gluon 
itself, because Z P (Q + ,Q 2 ) — 1, e q. (143), and the integral term in ( 166 ) vanishes. For K\ < Q 2 , the formfactor 
Z P (K + , K 2 ^) decreases (c.f. eqs. ( |162| ), fll65| )), and so the first term, which is the probablity that the gluon remains 
in its bare initial state, is suppressed by Z p , whereas the integral term, which is adjoint probability that the gluon 
contains a distribution of intrinsic gluons, increases with weight Z P (K + , K^) / Z P (K + , q\jz). Hence the evolution 
of the spectral density p describes the change of structure of the initially bare gluon state due to real and virtual 
emission and absorption of daughter gluons, corresponding to the generation of virtual Coulomb field coat and real 
brems-strahlung, respectively. 

Eq. (166) can be solved in closed form by using the following trick to effectively eliminate Z p . First, note that p 
satisfies the momentum sum rule Ifl9f , 



dq + q + p(q + ,q]_ 



(q\ql) = {K+f 



r 

JO 



dq + , 







,ql) = 0, 



(167) 



for any value of q\. Eq. (167) is nothing but a manifestation of lightcone-momentum conservation, meaning that the 
aggregat of q + momentum from intrinsic gluons must add up to the total K + of the gluon stat e com posed of those. 
This is a general property, which is immediately evident in the free- field case. Next, multiply (166) by q + / K + and 
integrate over q + from to K + , which yields on account of the sumrule (167), 



1 = Z p {K + ,K\) 



dq]_ a s {q] 



K 2 



2tt 



dzj{z) Z p 



K 



(168) 



which does not contain p. Next, multiply this formula with Z p 1 (K + , K]_) from the left, and then differentiate with 
respect to K\ by applying K\d/dK\: 



d 



Z p -\K\Kl)) p{K\Kl) + Z p -\K\Kl) {k\^p{K\K 



-Z» -\K\Kl)^p± f dzi{z) l V (^l, zK l 

Z7T Jq Z \ Z 



(169) 



Using (162), the derivative dZ p /dKj_ on the left hand side can be rewritten as 



Kl Mi ZP = - ZP - 1 {K+,Kl)\ asi ^- 



2tt 



dz j(z) 



(170) 



Substituing this into ( |166| ) and multiplying by Z p , one obtains a differential evolution equation a la DGLAP |B4 35 
that involves only p, but not Z p anymore FJ: 



K 



d 



■p{K+,K\) = 



a s {K 



1 dz 
— l\ z ) 



dKl ry " '"- L/ 2tt 
The explicit solution of this equation is well known |3qj3 



^Kl)- Z -p{K\Kl) 



(171) 



It should be noted that in obtaining (171), the fact that p 
to the very weak z-dependence of the the A^-argument of p pi 



K+ 



zKi 



K+ 



K\ J was used - a property that is due 
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p(K+,Kl) = Po (K+,Kl) + Pl (K+,Kl) exp 



where 



12tt v ±) 



exp 



^g(Kl) h(K+) 

117T 



p (K+,Kl) = 
Pl (K+,Kl) = 
g(K 2 ± ) = In 



(2tt) 4 
K+ 



- Q+) <5(K 2 - Q 2 ) 



ln(g 2 /A 2 ) 



ln(K 2 /A 2 ) 



47T 



1/4 



-3/4 



In 



(172) 



(173) 



F. Solving for the correlation function A^fJ 



Within the perturbative scheme ( |119| ) and (12C), the calculation of A cor — A^ r + ^Aj—j" is most conveniently split 



into two steps: 



1. The quantum contribution Aj^ r is evaluated to order g 2 from (107), i.e., the hard polarization tensor Tl cor [g 2 , A ] 



is to be calculated again in one- loop approximation with free- field internal propagators. The mean-field part 
polarization tensor, on the other hand, is set to zero in this first step: II =0. 

2. The mean-field induced correction SA^ in leading order gA is then added by calculating n c ° r [g; Ao]. The quantum 
part now is set to zero in this second step: n cor = (as it is already contained in AS? 1 " from Step 1). 



1. The contribution A% r 



Since to order g 2 only the radiative self-interaction contributes to to the hard propagator A^, and scattering 
processes that could alter the gluon trajectories are absent, the transport equation for the part simplifies to, 

0(9') ■ 



K-D r , A- 







(174) 



Therefore, with respect to the space-time variable r, (174) implies free-streaming behavior in the presence of the soft 



mean field, as implicitly contained in K = k — gA, that is, A^ r (r, K) = A c ^(r' — Vr ) with = K fl /K + and 



r' < r. Hence, it remains to consider eq. (|107|) with A 



Acor . 
^0 ■ 



a cor 

A [o] 



-A" f (a^ - 1 - n COT -[ 5 2 ,A ]) A« dv 



(175) 



The easiest way to obtain A^ is to use the formula (108) and simply convolute the total number density of gluons 



G = 1 + 2g with the spectral density p obtained in the preceeding subsection. To prove that the relation (108) is 



indeed consistent, one calculates instead Aj^ r from (175) directly. The procedure is fully analogous to the previous 



su bsect ion, except that, instead of IF et — IP , one needs to evaluate II > + II < . The resulting form of n cor [g 2 , Ao] 
in ftPnfy is 



IT 



2 



ip + iri (r,K) 

d 4 q 



""' '' (r,K) = - (i 

"" / 7^1 V 0l &(K,-q,-K + q) V Q d '^, v (r; K - q,q, -K) 

S ce' S M d xx '(q)d™'(K-q) {A>(r,q)A<(r,K-q) + A< (r, q) A> (r, K - q) } , (176) 
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where in the integral, the free-field forms Ag ^ v are given by (c.f. Appendix F): 

A %?(r,if) - Al(r,K)d^(K)S ab , A*(r,K) = {-2m) S(K 2 ) ( 0(±K+) + g (r, ±K)) . (177) 



Inserting into (176) the expressions (177), and observing that in ( |l07| ) n cor is sandwiched between A ret and A , 
i.e. appears only in the combination 



A^IT^A^ oc d^{K)d XT {q)d Tp {K-q)d pv {K) ex d^(K)U c r(K) 



(178) 



where n™ r = g pT H c p ° r , and d p \d\ u = d pvi one finds after a calculation analogous as in the preceding subsection the 
following result for the A^ part of the correlation function: 



A- (r,K+,K 2 x 



Tr 



dK'j 2 asiK'j 2 ) f 1 dz 



Kl K' 2 2ir 



J ^ 7 (z) [{p(l + 2 g )}(r,K + /z,K' ± 2 ) |{p(l + 2g)}(r,X+,if; 2 )] . (179) 



Comparison with (171) reveals that A^ r is indeed the convolution of the spectral density p with the total gluon 



density G = 1 + 2g, as advocated by eq. (|108j). Hence, 

dK' 2 a s { K ' 2 ) ' ' 



G n {r,K+,Ki) = G (r,Q+,0 ± ) 



k 2 , K' 2 



2tt 



dz 7(2:) 



K+ 

z 



l -G n (r,K+,K' 2 ) 



(180) 



where Gc^r, Q + , OjJ is the initial gluon density (154). In the limit z <C 1 the integral ( pL80|) can be approximately 
evaluated analytically p9| , pq| . This gives an estimate of the gluon multiplicity |4(],[4l| as a function of K\ at fixed 
space-time point r: 



r Q+ 

G m (r,K 2 x ) = i / dK+G m (r,K + ,K 2 ± ) 

J 



\\n{Kl/h?) 



exp 



Nc_ 
Hit 



'In 
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A 2 



(181) 



where Go(r, Q 2 ) is given by (154). It is evident that in the kinematic regime Q 2 > K\ > p 2 , the hard gluon multiplicity 
is characterized by a rapid growth as the gap between the initial scale Q 2 and K\ increases. 



2. The contribution 5A^' 
[A] 



The leading-order mean-field contribution ^Aj^ is now to be added to the result for A^ r , eq. (179). To do so, one 

■cor 



needs to evaluate II to order gA, using the free-field solutions (111)-(113) and set Tl cor — 0. The analogon of (175) 

(182) 



for SA'gT is 

[-4] 



6A m 



-A-'n [g,A ]A a ™ 



and II [g, Aq] can be read off from (104), giving the contribution 



9KxF X \r) 9fA- f (K) + g (r) A^(K) - A^(K) F Xv {r) 



ah 



gS ab d plJ {K) K x F M {r)d^A^ ah {K) . 



(183) 



32 



The second term on the left side cancels, because F X , A 



m— (o) An ~ F »9\v ~ ^ (nA-fiT« + K x n„)/K+, and ?i+ = -n = 
1, = 0, X - = 0. Notice that this is a specific feature of the employed lightcone representation, and does not hold 



(184) 



in a general non-covariant gauge. With fll83j ), the function SA^ satisfies the transport equation 



K-D r , SAffi 



g V K x F Xa d«A™ r 



To solve ( |184p , it is convenient to express 8A C ^ ^ in terms of a new function $ M = T a $^ 0, defined by 

6A<^Jr,K) = d, u {K) g**(r,K) d« A^(r,K) , 



[o] 



where AS? r is the solution (|179| ). In terms of the function <£ M , the transport equation (184) becomes now 



[o] 

The function $ M evidently satisfies 



[K-D r ,$^r,K)] = F^{r)K v 



if i $i 



, 



(185) 



(186) 



(187) 



i.e., $ is not an independent variable, but is expressable in terms of the transverse components and $ + is 
suppressed by K~ / K + and therefore may be set to zero. The interpretation of the function as was pointed out 
by Blaizot and Iancu @ , is the following: The component gtf?^ corresponds to the kinetic momentum K ^ = fc M — gA^ 
that is acquired by a gluon propagating in the presence of the soft field A^, or F^. The condition (187) reflects then 
the fact that the lightcone energy transferred by the soft field, namely g§>~ , equals the mechanical work done by the 
Lorentz force gV x ■ AK X = gVj_ ■ where V» = K^/K+ = (1, 0, V±) is the velocity. 

The transport equation (186) for $ M can be readily solved with the help of the retarded and advanced functions 

ret ret 
A adv J A adv 



A r ^(r-r',K)F^(r')K v - i 



A% dv (r-r',K) F^{r')K» 



(188) 



The free-field retarded and advanced functions admit the space-time representation 



A et (r - r', K) = -i 6( r - - r ~) S (r+ - r + - ( 



)) S 2 ( f± -r ± - |^(r -r ) 



A% dv (r-r',K) = +i 6{r ~ - r)s(r' + -r + - {r'~ -r - )) S 2 r[ - ?± - ^(r 



and therefore i (A r et - A$ dv ) (r -r',K) = 2 5 (r + - r'+ - (r~ - r'-j) 5 2 (V ± - r_[ - (K±/K+)(r 
tion into (p~8q) then yields, 



«!>,,( r.K) = 2 K" I <lr F,,,. ( r - -|- r 



2K V F^(r) 



(189) 
r ~~ ) ) . Inser- 

(190) 



Substituting this result into (185) and using the lightcone components of F M „, eq. (128), the result for the mean-field 
induced correction 6A?~ is, 

[A] 

6A^(r,K) = T r \d-J(K)SA^; n Jr,K)] = -gF ±+ (r) ( K ± 



[A] nu 



dK+ 



- K+ 



dKj 



A^(r,K) 



(191) 



where _L denotes the transverse vector components i — 1,2, and &± — ^( < i )1 + "i" 2 ), K± — -^(K 1 + K 2 ). 
With jl90|), the addition SG^ (r, K+,Kj) to the gluon density G m (r, K+, Kj_) of eq. (p9) is, 



K+ _ 



d 



8G^(r,K+,Kl) = -2.9 — F ±+ (r) (X* Gpj,(r, K+ ,K\ ) J + O (Xi/(A'+) 2 ) , 
where the explicit form of derivative term in brackets can be easily read off the right-hand-side of (|l80| ) 



(192) 
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G. Expansion in space-time of the hard gluon ensemble 



As argued before in (|17j), the evolution of the hard gluon density G, described by A cor , can in the present context 
be viewed as a purely multiplicative casade of gluon emissions , since to order g 2 and due to the quasi-collinear 
motion of the gluons, statistical scatterings between them do not contribute. Therefore the space-time development 
of G(r, K + , Kj_) with respect to r = (r~, r + , is of 'free-streaming nature'. That is, the expansion with time 
of the ensemble of gluons as-a-whole proceeds through a deterministic diffusion in momentum and space-time, as 
qualitatively scetched in Fig. la. 

To quantify this heuristic picture, one needs to invoke the uncertainty principle to relate the development in 
space and time to the evolution in momentum space, i.e. with respect to K + and K\. Specifically, what is the 
characteristic time in the chosen Lorentz fame, that it takes to build up the density G(r, K + , Kj_) from the 
initial form Go (70, Q + , 0±) at time r$ = Q. Viewing the gluon evolution as a cascade of successive branchings 
K n -i — > K n + K' n , where n labels the generation in the cascade tree, the life-time of gluon K n -\ is given by the time 
span Ar~ that it takes to emit and form the daughters K n and K' n as individual off-spring, that is, by the formation 
time 



Arl 



Kin 



K' 2 

^±n 



K- 



K 2 

^±n 



T n In 



(193) 



with K+ = K+_ x - K+, Kj_ n = -K' ±n , and K ±n = \J K\ n . Here t„ = 1/K± n and j n = K+_JKx n play the role of 
the proper time and the Lorentz gamma-factor, respectively , in agreement with the uncertainty principle. Similarly, 
the average longitudinal and transverse distances travelled by the gluons K n and K' n during the time span Ar~ are 



Ar 4 

1 — T-i 



Arj_ ri 



V + + V 

r n 1 y r 



+ ' 



v ±r , 



Arl 



K 2 

^±n 



where V^ 1 = /K + and K± <C K + is assumed as before. The average total time ( 
cascade generation with mean gluon momentum K + and K^, and the associated spatial spread (r + ) , (r ±), of 



(194) 
(195) 

elapsed up to the n-th 



the diffusing gluon ensemble, is then o btain e d by weighting the evolution of the gluon den sity G, eq. ( 180 ), with 
Ar(K+,Kj) ee (At- - , Ar+, Ar ± ) from (^93|)-(^95|). Taking the real emission part of eq. ( |l80| ), differentiating it with 
respect to K±, convoluting it with the weight Ar(K + , Kj_), integrating over all possible branchings, and normalizing 
it to the density G(r, K + , Kj_) itself, the desired average is 



r{K + ,K\) 



G(r,K+,Kl) 



q2 dK 2 a s {K' 2 ) 



K 1 , 



K' 2 



2tt 

- 1 dz 



dK- 



G(r,K+ ,K'_?) Ar(K+',K 



K + 



' 2\ 



J ^ 7 WG(r,X+/2,X; : 



(196) 



This complicated formula can be approximately evaluated from the known behaviour of G, as has been worked out in 
detail in fl2fl. For K + <C Q + the result is, up to powers of \n(Q + /K + ), the following estimate: 



( r-(K+,Kl))~( 



-(K + ,Kl) ) = ^T(K^ 



r ± (K+,K 2 ± ) ) = — T(K+,K 



where 



(197) 



T(K+,Kl) = Cl (Kl) exp 



(198) 



with ci , C2 > very slowly varying functions of K\ . This estimate shows that those gluons which are emitted either 
with large K^, or with small K + /Q + , appear the earliest in time r~ and contribute the quickest to the diffusion in 
r+, rj_. 



34 



H. Constructing the hard current j M and the induced soft field 



The final task of the solution scheme Sec. Ill C is to solve for the soft field A^, or F M „, which is induced by the 
color current j M , being generated by the aggregat of initial plus emitted hard gluons from the evolution of the gluon 
density (18C). In the equation of motion for F^„, recall (105), 



(r) = 



the current on the right-hand-side is determined by the hard gluon correlation function A™'' — Af-, r 



and therefore by the gluon density G = G 
%{r) = T a %{r) = -g 



[0] 1" SG [A]' 

Tr 



as obtained in the previous subsection, 
T a (k^A^'^K) - K v A™ rv {r,K) 



(199) 



[A] nv' 



(200) 



The first point to be made here is that, for the lightcone gauge condition A + = 0, the gauge-field tensor F ^ 
has only the non- vanishing components (126), and if one requires in addition A~ = 0, then in (|l99| ), D X — 



8 j_ 5^+ D± Fj_ + . The second observation is that the left-hand-side of (200) is essentially the density G of hard gluons 



weighted with their momentum . Because the gluons evolve with the velocity = K'- l /K + along the lightcone, at 
a given lightcone time r~ and corresponding coordinate r + = + V + r~ = r~ , these gluons appear as an extremely 
thin Lorentz-contracted sheet, but are spread out in transverse direction r± over a disc with radius ~ 1/(K±). As a 

consequence, the gluon current j M = Q + , j~ , j ±) has only a component in +- direction, j M = j + [^J,^3|. Denoting 
as before the two transverse vector components i = 1, 2 by _L with summation convention a±b± = J2i=i i a ibi, eq. 
( |199| ) now becomes, 



D±, F +± 



(r) 



^- 5 r bc Al(r), Fl x (r 



J+(r) 



where is the color-charge density at r = r + , 

f + {r) = gT a J(r ± )5(r--r+) 



where 



J{ri_) = 2ir 



and, A cor = d~l (A^T + SMS ) 

V v v [0] ytv [A] iiv' 



(2i:f2K- ,,, 

using pa, pd, 



dKi Tr 



(201) 



(202) 



(203) 



A co '(r,X+,ifi) =2(1- 2g—F ±+ (r)K 



OKI 



(204) 



Eqs. ( p03| ) and (204) follow from the fact that, on the left-hand-side of ( |200| ) , the correlation function obeys the 
transversality condition K^A^J = K^d^A cor = and because K»A c ° r v = K^d^A cor = 2K^A cor . Notice that in 
(203) the limits of the integra tion over K + and K\ correspond to the average time r~ and spatial extent r± of the 



gluon system, as estimated in (197) above, and hence, J accounts for the total gluon multiplicity accumulated by the 
evolution between Q and /i. 



Integrating both sides of (201) over r + , r , and using (|203D, (204), gives 



J(r±) = dr+ dr- 

Ja Jo 



- ig[A ± ,d ± F ±+ ]){r+,r-,r ± ) 



(205) 



An approximate method to determine the soft field from (205), is to adopt the approach of Kovchegov M3j, who recently 
calculated the lightcone gauge field induced by an ultra-relativistic current of quarks with a uniform momentum 
distribution, using the known form of the lightcone gauge potential of a single color charge |44| . Applying his concept 
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to the pre sent case of gluons with a non-uniform distribution G(r, K), the first step is to write the color-charge density 
J of ( 203 ) in a 'discretized version' as a superposition of M individual gluon charges, 



j a +(r) = gTr 



AT 



rf)S(r- - rr)S 2 (r 1 _ - rli) J(f ±i ) 



(206) 



where AT is the total number of gluons at a given r + — r , Af(r) — tt (2*y'2K+ J/j? dK± G(r,K). Now the 
approximate solution to ( p05[ ) for the lightcone gauge-potential at space-time point r is obtained, by the superpo- 
sition of contributions that are induced by the hard gluons at points r j . Following Kovchegov Ji3| in detail, the result 
is that has only non-vanishing transverse components, 



A + (r) 



A ± (r) 



A~{r) = 

2ng 0(r + -rt)< 



In 



r_L-r_L, 
\r± - r ±i 



Jlfxi) Tr[Tf S{r)T a S- l {r)] 



where 



Sir) 



JJexp 



2n g 2 T a T? 6{r+ - r+) 6(r~ - rT") In 



r±-rj_, 



(207) 



(208) 



It is important to note that eq. (207) is only an approximate solution of (205) for the induced soft field. It is an 

estimate of the classical equation of motion for the soft mean field A± that is generated by th e collective motion 
of a given configuration of hard gluons with a distribution G(r, K). In other words, eq. (207) is the non-abelian 
Weizsacker- Williams field due to the hard gluons. 



I. Resume 



Let me summarize the input and results of the preceding sample calculation for the evolution of a high-energy gluon 
beam along the lightcone. On the basis of the calculation scheme of Sec. Ill C, the logic of application proceeded in 
the following steps: 



Choice of lightcone gauge with gauge vector n along the gluon beam direction K + /K and gauge constraint A^ 



0. 



Specification of the initial bare gluon ensemble at time r~ — with a momentum distribution of equal momenta 
K + = Q + , K\_ = 0, and a spatial distribution being uniform rj_ < R in the transverse plane, but a delta- function 
sheet in longitudinal beam direction at r + = 0. 



3. Calculation of the retarded and advanced functions Aadv and the associated sp ectr al dens ity p to order g 2 from 
the initial values of the hard gluon propagators. The result is stated by eqs. (172), (173). 



4. Evaluation of the quantum part Aj^ r of the correlation function, involving the result for Aadv of point 3. The 



solution for A22 r and the corresponding gluon phase-space-density is given by (180) and (181), respectively. 



[0] 

5. Evaluation of the mean-field part SA^, involving the solution for A? 

[A] 



correction to the gluon density SGrj, is given by (191) and (192), respectively. 



of point 4. The result for SAf— * and the 

[A] 



6. Construction of the hard gluon current j from the solution A cor = A 



[o] 



^ApT of points 4. and 5., with exploit 



form given by the formulae ( ^02[ )-(204). Approximate evaluation of the soft mean field A from the classical 
Yang-Mills equation (201) with resulting Weizsacker- Williams form ( 207 ). 



With this procedure, the original master equations (102)- (105) are solved in first iteration to order g 2 (l + gA). One 
could in principle now repeat this cycle, with the first order solutions replacing the zcroth-order forms as input. This, 
however, is another story. The story of this paper ends here. Cheers. 
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APPENDIX A: CONVENTIONS AND NOTATION 



Throughout the paper pure SU(3) C Yang-Mills theory for N c = 3 colors is considered, in the absence of quark 
degrees of freedom, with the gauge field tensor 

K» = d l A l - d l A l + 9f abc AlAl , (Al) 
and the classical Yang-Mills Lagrangian 

C YM (x) = \j%,{x)^' a {x) = - \ { {dlAlf - (d*Al) (d^ a )} (x) 

+ 9 fabc { {d?Al) A^ b A v ' c } (x) + g 2 f<*°f°M { A\A C V A^ b ' A v > c ' } (a;) (A2) 

Because only gluonic degrees of freedom are considered, only the fundamental representation of color space is relevant, 
with the color indices a, b, . . . running from 1 to JV C . The generators of the SU(3) color group are the traceless hermitian 
matrices T a with the structure constants f abc , as matrix elements, satisfying 

Tr (T a , T b ) = N c S ab , [T a ,T b ] = +i f abc T c , -i f abc = (T a ) bc . (A3) 

In compact notation, 

A/j — T A^ , -Fftv — T J- d^Au d^A^ ig \A^, A v \ ,D U ] > (A4) 

where = d/dx^ acting on the space-time argument — (x°,x). The covariant derivative, denoted by D^, is 
defined as 

£>„(*) = 31 - igT a Al(x) =8%- igA^x) , (A5) 

and its adjoint is D^(y) = + igA^(y). In components, using (|A3|), 

Df{x) = 5 ab d? - gf abc A°(x), (A6) 

with the color coupling strength g being related to the strong coupling a s = g 2 / (47r). In general, for any color matrix 
O with matrix elements O a b{x), the action of the covariant derivative is 

[D^O{x)\ = d*0(x) - ig [A„(x),0(x)} , (A7) 

and in particular, the covariant derivative of the field strength tensor reads [D^ ,F V \\ — d^T v \ — ig [A^^T^x]. 

The convention for placing indices and labels are such that color indices a,b, . . . are always written as superscripts, 
whereas all other labels may be subscripts or superscripts. In particular, the Lorentz vector indices v, . . . may 
be raised or lowered according to the Minkowski metric g^ — diag(l, — 1, — 1, — 1), and the usual convention for 
summation over repeated indices is understood. Finally, some shorthand notations are employed, namely 

A ■ B = A^ V B V , K ■ (A B) = A»B V (A8) 

Ao B = J d i xA{x)-B(x), Ko(AB) = J d A xd A y K(x,y) ■ (A(x) B(y)j , (A9) 

where the label P under the integral sign refers to the integration of the time components x° (y°) along a closed path 
in the complex time plane. 



APPENDIX B: BASICS OF THE CLOSED-TIME-PATH FORMALISM 

1. The in-in amplitude Zp 

The key problem in this paper is to describe the dynamical development of a multi-particle system (here gluons), 
that evolves from an initially prepared quantum state, e.g., produced by a high-energy particle collision. There is 
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a crucial difference between the evolution of the system in in vacuum (which means, free space in the absence of 
surrounding matter), and in medium (which could be either an external matter distribution, or an internal particle 
density induced by the gluons themselves) . As illustrated by Fig. 2 in the introduction, this difference arises from the 
interactions, and hence, non-trivial statistical correlations between the gluons and the particles of the environment. 

In the case of vacuum, the usual quantum field theory desrcibes the time evolution of the system by the vacuum- 
vacuum transition amplitude, also called in-out amplitude (see Fig. 2a, left panel). That is, one starts at time to in 
the remote past with appropriate asymptotic m-state and evolves it to too in the asymptotic future, by means of the 
time evolution operator U (too, to)- Multiplication with the hermitian conjugate counterpart, which corresponds to a 
backward evolution from too to to under the action of (to, too)- The resulting in-out amplitude may be interpreted as 
the sum over all n-point Green functions for space-time points along a path in the complex t plane, exclusively on the 
upper (lower) branch for the forward (backward) evolution. In vacuum there is no correlation between the two time 
branches, and so, for instance, the 2-point Greenfunctions are the usual time-ordered Feynman A F (anti-time-ordered 
A F ) propagator (sec Fig. 2a, right panel). Because W (to, too) = U(too,to), one has A F (ti,t 2 ) = —A F (t 2 ,ti). 

In the case of a medium, the above concept fails, because of the a-priori-presence of medium particles described 
by the density matrix p(to)- Instead one has to construct a generalized transition amplitude, called in-in amplitude, 
which accounts for the non-trivial initial state at to embodied in the density matrix p(to), and evolves the system 
in the presence of the medium from to to too in the future, by means of the time evolution operator U (to, too) (see 
Fig. 2b, left panel). Because now W (to, too) 7^ U(too,to)p(to), forward and backward contributions are not merely 
conjugate to each other, but interfere, giving rise to statistical correlations between upper and lower time branch of 
the contour in the t- plane. As a consequence the space of Green functions is enlarged by non-time-ordered correlation 
functions. For example, the 2-point functions are now A F , A F plus the new functions A< and A> (see Fig. 2b, right 
panel) . 

The fundamental quantity of interest is the in-in amplitude Zp for the evolution of the initial quantum state \in) 
forward in time into the remote future, starting from a specified initial state that could be either the vacuum or a 
medium. Within the CTP formalism the amplitude Zp can be evaluated by time integration over the closed-time-path 
P in the complex t-plane. As illustrated in Fig. 6, this closed contour extends from t = to to t = too m the remote 
future along the positive (+) branch and back to t — to along the negative (— ) branch, where any point on the + 
branch is understood at an earlier instant than any point on the — branch: 

Z P [J,p] = Z P [J+,J-,p] = Tr{uy(t ,t)U(t,to) J+ p(to)} , (Bl) 

where the p(to) is the initial state density matrix, and U and W are the time evolution operator and its adjoint, 

Uj+{t,t ) =Texp|-i jf dt'd 3 x'J + (x')-A + (x') | 

U ] J+ (t ,t) =T f expj+i jf dt'd 3 x' J~(x')- A~(x') J , 



(B2) 



with T (Tt) denoting the time (anti-time) ordering operator. Note that J + (J_) is the source along the positive 
(negative) branch of the closed-time-path of Fig. 6, and in general J + ^ J_, so that Zp depends on two different 
sources. If these are set equal, one has Zp(J, J, p) = Tip, which is equal to unity in the absence of initial correlations, 
being a statement of unitarity. 

Zp contains the full information about the development of the initial state via the creation, interaction, and 
destruction of quanta, through the agency of the sources: the quanta are initially created (e.g., by particle collision), 
they evolve by further creation and annihilation (real and virtual emission/ absorption as well as scattering), and are 
finally destroyed (e.g., by detection in a calorimeter). Both the act of initial creation and final destruction represent 
the external sources J in the sense of a probing aparatus, whereas the intermediate dynamics is governed by the 
underlying quantum theory. Hence, in order to describe the time evolution of the initially prepared quantum system, 
to the final detected state, the knowledge of Zp allows to extract objectively the self-contained development of the 
system, when the external influence removed (i.e. the sources are switched off). 
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FIG. 6. The close-time-path in the complex t-plane for the evolution of operator expectation values in an arbitrary initial 
state. Any point on the forward, positive branch to — > too is understood at an earlier instant than any point on the backward, 
negative branch too — > to. 



The functional Zp can be represented as a path integral by imposing boundary conditions in terms of complete 
sets of eigenstates of the gauge fields A^ at initial time t = to, 



A(t ,x)\A + (t )) = A + (x)\A+(t )) , 
and in the remote future at t = t^ , 

•4'(too,£)|.A'(ioo)) = A'ffilA'itn)) 



A(t ,x)\A-(t )) = A-(x)\A-{t )} 



(B3) 



(B4) 



Then, making use of the completeness of the eigenstates, one obtains from (Bl) the following functional integral 
representation for Zp, 



Z P {J+,J-,p}= J VA+VA-VA 1 {A-(t )\uy{t ,t oo )\A{t oc )) 

x (A(t oo )\U J+ (t oo ,t Q )\A + (t )) x (A + (t )\p\A + (t Q )) 



(B5) 



The first two amplitudes are the transition amplitudes in the presence of and J~ , whereas the density matrix 
element incorporates the initial state correlations at to at the endpoints of the closed-time path P. Hence, one obtains 
the path integral representation for Zp in analogy to usual field theory |2Q,|2l| 



Zp[J+,J-,p} = J VA+VA- exp[i(l[A + ] + J + o A + ^j - i(l*[A~] + J~ o A~) M[p] , 

where 

M(p) = (A + (t Q )\p\A-(t )) . 
The generalized classical action I [A] accounts for all four field orderings on the closed-time path P: 



I[A] = I[A + ] -I*[A~] = I^[u a0 AlA v ] + ^[gv^idMA^Al] + I^[g 2 w^AZA^AlA 



5 ] 



(B6) 



(B7) 



(B8) 



where the correspondance with the terms with the ones of i 
where a, /3,7, 8 = +, — , 



is obvious (the color indices are suppressed here), and 



laP = u afi = diag(l,-l) 



"q/3 7 = Sapup^ , w a p lS = sign(a) 8 a0 8 l3l 8 lS , 



(B9) 



with the usual summation c onve ntion over repeated greek indices a, /?,... Eqs. (B6)-(B8) represent the detailed 
version of the compact form (CI) used in Sec. 2 as the starting point, except for the Fadeev- Popov determinant and 
the gauge fixing constraint, which is omitted here, and which is addressed in Appendix C. 
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2. The density matrix and the initial state 



Turning to the properties of the initial state incorporated in the functional M(p) with the density matrix p(to), we 
denote by to the initial point of time from which on the evolution of the multi-gluon state is followed, and assume 
that all the dynamics prior to to is contained in the form of the initial state 



\in) = \A(t )) = II 



(BIO) 



The initial state at t — to can be constructed by expanding the gauge field operator A in the Heisenberg representation 
in terms of a Fock basis of non-interacting single-gluon states, the m-basis, 



d 3 k 



(Bll) 



so that a particular Fock state is given by (supressing color and Lorentz indices) 



,n 



(oo) \ _ 



VnW! 



(ki-, Si) 



(B12) 



Here & (c) is the creation (annihilation) operator for a single-gluon state with definite 4-momentum k and spin s, 
satisfying # | 0) = 1 (e| 0) = 0). and the nfi are the occupation numbers of the different gluon states, nW = 
{n^'\c(ki, Si)c^(ki, s<) | nW ). Finally, |0) denotes the vacuum state or a ground state different from vacuum (e.g. a 
hadron). Thus, a general multi-gluon state \<j>) at time to is given by a superposition of such states, 



\A(t ) 



,(«=)■ 



(B13) 



A», a ri(0 



with real-valued coefficients C. Alternatively, the initial state of the system at to can be characterized by the density 
matrix, 



p(t ) = \A(t ))(A(t )\ 



(Mi 



(B14) 



For instance, the case of empty vacuum corresponds to a diagonal density matrix p(to) = |0)(0| with (po)^ oc Sij, 
whereas a general density matrix that describes any form of a single-particle density distribution at t Q is 



p(to) = N exp 



E 



d 3 k 



(2?r) 3 2fc ( 



■e(k )F(t ,x,k) c*Hk,s)c«(k,s) 



(B15) 



where Q denotes the hypersurface of the initial values and F is a c- number function related to the single-particle phase- 
space density of gluons around x + dx with four-momentum within fc M +dk°dk, and iV a normalization factor. The form 
(B15) describes a large class of intreresting non-equilibrium systems fl9| , and contains as a special case the thermal 

equilibrium distribution, namely when t — > —i/T and F{to, x, k) — » fco <5(fc 2 ) T _1 , so that /5(Tq) — > TV exp —Hym/T 



3. Perturbation theory and Feynman rules 



The convenient feature of the CTP formalism is that it is formally completely analogous to standard quantum field 
theory, except for the fact that the fields have contributio ns fro m both time branches. In particular, one obtains as 
in usual field theory, from the path-integral representation (Cll) the n-point Green functions G^ n '{x\, . . . , x n ), which 
however now include all correlations between points on either positive and negative time branches, 



^aia 2 ... a„ (-"li ■ ■ ■ ) x n) — 



1 



Z P [0] idJC^ 



-Z P [K] 



± 



(B16) 



K=0 



depending on whether the space-time points Xi lie on the + or — time branch. One can then construct a perturbative 
expansion of the non-equilibrium Green functions in terms of modified Feynman rules (as compared to standard field 
theory) ||M||. 
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(i) All local 1-point functions Ga\x), such as the gauge-field or the color current, are 'vectors' with 2 components, 

A(x) EE W = ( J-) ( B17 ) 

Similarly, all 2-point functions G^p(x, y), such as the as the gluon propagator iA^ v and the polarization tensor 
11^, arc 2x2 matrices with components, 

A(xi,a:2) = ^ A -+ A — J n(x u x 2 ) = ^ n _+ n — J • (B18) 

Explicitcly, the components of the propagator are 

A^(x,y) =A++(x,y) = -i ( T .4+0*0 A+(y) ) A<„(x,y) = A+"(a;,y) = -i ( A+(y) A~(x) ) 

A>„(x,y) =A-+(x,y) = -i { A~(x) At (y) ) tf v (x,y) = A~(a:,y) = -i ( T^(a:) ^"(») ) , (B19) 



where A f is the usual time-ordered Feynman propagator, A F is the corresponding anti-time-ordered propagator, 
and A> (A < ) is the unordered correlation function for x > yo {xo < Uo)- In compact notation, 

A^(x,y) = -i(T P A(x)A(y)(y)) , (B20) 

where the generalized time-ordering operator Tp is defined as 

T P A(x)B(y) := P (x o , y ) A(x)B(y) + 9 P (y , x ) B(y)A(x) , (B21) 



with the 6*p-function defined as 

a , s J 1 if .to succeeds yo on the contour P , 
P (£0,2/0 J — if ,x precedes j/o on the contour P ' ^ ' 

Higher order products A{x)B(y)C(z) . . . arc ordered analogously. Finally, the generalized i5p-function on the 
closed-time path P is defined as: 

+5 4 (x — y) if xq and yo from positive branch 
5 P (x,y) := ^ —5 i (x — y) if x and yo from negative branch . (B23) 
otherwise 



(ii) The number of elementary vertices is doubled, because each propagator line of a Feynman diagram can be either 
of the four components of the Green functions. The interaction vertices in which all the fields are on the + 
branch are the usual ones, while the vertices in which the fields are on the — branch have the opposite sign. 
On the other hand, combinatoric factors, rules for loop integrals, etc., remain the same as in usual field theory. 
Specifically, the 3-gluon and 4-gluon vertices, 

G%( Xl ,x 2 ,x 3 ) ee J p d 4 xG^ a (xi,x)Gf) (x 2 ,x)G^(x 3 ,x) T% 7 ,(x) 

G { ^ 5 {x 1 ,x 2 ,x 3 ,x 4 ) = J^d i xG^)jx 1 ,x)G ( ^{x 2 ,x)G^) 1 {x 3 ,x)G^ ) s {x4,x) T^) fj , ys ,{x) , 

with r( 3 )(:r) and r( 4 )(a;) denoting the elementary, amputated vertices (with the external legs removed), have, 
for fixed a, (3, 7, S two components. For instance, as in Fig. 7, for the external points on the +-branch, 

for a = P = 7 = + : T% y (x) = (vf ++ , ifL) 
i ova = (3 = 1 = S = +: T%, ( , s ,(x) = (v^ +++ , T^___) 



42 



(P\xp x 2 ) = i A (xj , x 2 ): • • •- 



+ + + 



G (xj,x 2 ,x 3 ): "+ 





(j \xj, x 2 , x$, x 4 ): An 



+ + 



FIG. 7. Example for the appearance of additional contributions to the n-point functions G*-™' for the propagator G^ 2 \ the 
3- vertex G^ 3 ', and the 4- vertex G' 4 ' In usual quantum field theory referring to free space or 'vacuum', only the +-graphs are 
non-zero. In the CTP formalism, accounting for the presence of surrounding matter or 'medium', new diagrams arise that 
correspond to statistical correlations between the field living on the + and — time branches of Fig. 6. 

APPENDIX C: THE IN-IN AMPLITUDE Z P FOR QCD IN NON-COVARIANT GAUGES 

For the case of QCD, path-integral representation of the in-in amplitude Zp is obtained along the lines of Appendix 



Bl, except that one has to extend the generic formula (B6) to account for elimina ting the spurious gauge degrees of 



freedom by the usual Fadeev-Popov procedure ||45| . The gauge theory version of (BE) for the class of non-covariant 
gauges (||) reads therefore: 

Z P [J,p] = N J VAdetF 5(f[A}) exp |i ( I [A, J]j } M(p) , (CI) 

where A = (A + ,A~) and J = (J + ,J~) have two components, living on the + and — time branches of Fig. 2. 
Physical expectation values are defined as functional averages over Tp-ordered products of n field operators (n > 1), 
weighted by Zp, 

(0 1 (x 1 )...O n (x n )) P = p j VA detF 5 (f[A\) exp [A, J})} M{p) Tp [Oi(xx) . . . O n (x n )] . (C2) 
The structure of the functional Zp in ( |Clj ) and ( |C2[ ) is the following: 

(i) The functional integral (with normalization Af) is over all gauge field configurations with measure T>A = Yin a ^-A^, 
subject to the condition of gauge fixing, here for the class of non-covariant gauges defined by 

f a [A] := n-A a (x) - B a (x) (n»A«(x)) =0, n" = -= , (C3) 

V\ n I 

where n M is a constant 4- vector, being either space-like (n 2 < 0), time-like (n 2 > 0), or light-like (n 2 — 0). With 
this choice of gauge class the local gauge constraint on the fields A a ^{x) in the path-integral ( ]Cl| ) becomes, 

detJ 7 5 (n • A a - B a ) = const x exp {i I GF [n ■ A}} (C4) 
I GF [n-A] = ~ [ d A x [n-A a (x)f , (C5) 



p 



43 



where det^F is the Fadeev-Popov determinant (which in the case of the non-covariant gauges turns out to be a 
constant factor, c.f. Appendix C), and where S(n-A) = Yl S(n-A a ). The right side translates this constraint into 
a the gauge fixing functional Iqf- The particular choice of the vector n M and of the real- valued parameter a is 
dictated by the physics or computational convenience, and distinguishes further within the class of non-covariant 
gauges g 

homogenous axial gauge : n 2 < a = 
inhomogenous axial gauge : n 2 < a = 1 
temporal axial gauge : n 2 > a = 

lightcone gauge : n 2 = a = . (C6) 

(ii) The exponential / is the effective classical action with respect to both the + and the — time contour, / JT'] = 

/ [A + , J + ] — I* [A~ ,J~], including the usual Yang-Mills action Iyu = J d A xLyu, plus the source J coupled 
to the gauge field A: 

I [A, J] = -\J^xT; v {x)T^ a (x) + J^xJ«(x)A^ a (x) 

= Iym[A] + Jo A. (C7) 

(iii) The form of the initial state at t = to as described by the density matrix p (an example is given in Appendix B2, 



eq. (B14)) is embodied in the function M.{p) which is the density-matrix element of the gauge fields at initial 
time to, 

M(p) = (A+(t )\p\A-(t )) = exp{iK[A]) , (C8) 

where A refers to the + and — time branch at to, respectively (c.f. Fig. 2). The functional JC may be expanded 
in a series of non-local kernels corresponding to multi-point correlations concentrated at t = to, 

K[A]=K.W + f d A xK^ l \ a {x) A"' a (x) + I f d^ylC^fix^) A"' a (x)A u ' b (y) ... 
J p 2 J P 

= E ^ W ^(l)^)-iW) ■ (C9) 

ri=0 

Clearly, the sequence of kernels K.^ contains as much information as the original density matrix. In the special 
case that p is diagonal, the kernels JC^ = for all n, and the usual 'vacuum field theory' is recovered. 



The path-integral representation (CI) can be rewritten in a more convenient form: First, the gauge-fixing functional 
3f[^'^4] is implemented, usi 
4(p). Third, K, (t)) is absorbt 
J in the 1-point kernel K.^: 



lGF[n-A] is implemented, using (C5). Second, the series representation (C9) is inserted into the initial state functional 
M.{p). Third, IC^ ' is absorbed in the overall normalization N of Zp (henceforth set to unity), and the external source 



/C (0) := i hxM , JC {1) := /C (1) + J. (CIO) 

Then ( |Cl| ) becomes, 

Z P [J,p] => Z P [K] = J VA exp{i(l[A,tcf)} , (Cll) 



where, instead of (|C7|) , 



§ § The analogy with the class of covariant gauges defined by / a [-4] := d x ■ A a — B a , instead of (|C3|), is evident: in place of 
(35), it results in the familiar gauge- fixing functional exp • -i/2a J p d 4 x(d-A a ) 2 )}, where a = 1 gives the Feynrnan gauge 



and a = the Landau gauge. 
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I[A,K] = I Y m[A] + I G F[n-A] + K {1) oA + ^JC i2) o(AA) + i/C (3) o (AAA) 



(C12) 



The objects of physical relevance are the n-point Green functions G*-™- 1 , defined as the coefficients in a functional 
expansion of Zp, 



oo . n p n 

Z P [K] = Z P [0] Hd 4 Xi G^( Xl ,...,x n ) JCW( Xl ,...,x n ) 

m — i '"'J ■ — i 



(C13) 



that is, the are functional averages in the sense of (|C2]), 



Z P [0] i6Kt n ) 



Z P [K] 



(C14) 



K=0 



The practical evaluation of Zp amounts therefore to calculating the G*-™- 1 in the ex pansi on (|C13| ) up to the order of 
desired accuracy. For instance, the 1-, 2-, and 3-point Green functions according to (C14) are 



G^;(x) = (A«(x))p 

G^t(x,y) = {A a Jx)At(y)) P = 



i8K, {1 l a {x) 



<(y))p + < A a Jx) ) P { Al{y) )j 



+ f ma, : ( < Ab M M A ^ Z ) )p) + < A » M Ab M )p{ Al(z) ) P . (C15) 

ioK 'fj, (x) v ' 



Higher order Green functions are generated similarly from ((C14) 



APPENDIX D: NON-COVARIANT GAUGES AND THE ABSENCE OF GHOSTS 

In this appendix the standard procedure of gauge field quantization is applied to the class of non-covariant gauges 



(C6), and it is shown that ghost degrees of freedom are indeed absent, reducing the general non-linear dynamics in of 
QCD essentially to a linear QED type dynamics. For an excellent review and bibliography, see Ref. p7| . Recall that 
under local gauge transformations 



g[6 a ] = exp(-i6» a (x)T a ) 



the gauge fields transform as 



AO. 



A 9)a = g[6 a ] At g- 1 ^] , 



(Dl) 



(D2) 



implying that T^T^ — T^J a T$ a , that is the gauge invariance of the Yang-Mills action Iym[A]. However, the 



source term J o A in the generating functional Zp of ( |Cl| ) is not gauge invariant under the transformations ( pl| ). 
Consequently, the naive functional 



4™ 3 = j VA CX P { l ( Iym [A]+J°A)} x M(p) (D3) 



is also not a gauge invariant quantity. As is well known, this can be remedied by applying the formal Fadeev-Popov 
p5[ procedure and integrate in the path-integral Zp over all possible gauge transformations g(0 a ) subject to the linear 
subsidiary condition 



4> a [A^} = n»A^ a (x) - /?» = 



(D4) 
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with normalized space-like vector n M and f3 a (x) an arbitrary weight function. The Fadeev-Popov trick to implement 
the constraint (D4) in the non-invariant functional z^ ame ) by multiplying with 



i = fve n^(^°[4 e) ]) dctF > 



where the determinant is the Jacobian for the change of variables 4> a — ► 6 a 



( detJ^ 



det 



ve l[sU a [A[ 



(D5) 



(D6) 



Following this procedure one arrives at 



7 (naive) 



Z P = J VA ActT JT 5 {<t> a [A^}) exp {i (l Y M [A] + J o A)} x M(/3) , (D7) 



which is now a gauge invariant expression due to the proper account of the subsidiary condition (D4) that guarantees 
the correct transformation properties of the gaug e fie lds in the presence of the sources J . 

To obtain t he f inal form of Zp as quoted in (CI), one integrates functionally over the arbitrary funtions (3 a (x) 
introduced in (D4), by choosing, e.g., a Gaussian weight functional 



w\j3 a ] = expj-^- J p d 4 x [f3 a (x)} 2 } , 



(D8) 



with the real valued parameter < a < 1, upon which the Fadeev-Popov determinant det J 7 can be rewritten in a 
more suitable way: 



detT 



exp 



± f d^x [/3 Q (x)] 2 l 6(n»AW a (x) - p a (x)) 



(D9) 



In order to calculate the determinant, it is sufficient to integrate over 9 a in a small vicinity where the argument of 
the 5- function passes through zero at given A^> a and [3 a . For infinitesimal gauge transformations 



g[o a ] 



the gauge fields transform as 



A a 



so that one obtains 

d(n^Af a (a 



A% + 5 A" 



a ] = 1 - iO a (x)T a , 



SA a = fabc 9 b A c _ lgx a 

g 



(D10) 
(Dll) 



/T(x)) = s(n"A^ a (x) + f abc 6 b n^A [ * )c - ^n^d*9 a - (3 a ^j 



1 



8 [ f abc 9 b f3 c - - n^d* 9 a 



(D12) 



because n^A^ ] a = (3 a . This latter expression is evidently independent of A®. Therefore, when substituted into ([OKI), 
det.F is also explicitely independent of the gauge fields, and hence can be pulled out of the path-integral Zp and 
absorbed in an (irrelevant) normalization, which may be set equal to unity. The final result is then: 



where, from (D9), 



Z P [J,p] = J VA exp{i(/y M [.A] + I G F[n- A] + J o A)} x M{p) 
I GF [n-A] = cx p|-^ J d i x[n-A a {x)] 2 ^ . 



(D13) 



(D14) 



In conclusion: the property of gauge field independence of the Fadeev-Popov determinant proves that there are indeed 
no ghost fields coupling to the gluon fields, hence the formulation is ghost-freel 
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APPENDIX E: THE TRUNCATED EFFECTIVE ACTION I> A, A 

The generating functional for the connected Green functions, denoted by Q^ n \ is defined as usual: 

W P [K] = -i In Z P [JC] . (El) 



From Wp one obtains the connected Green functions by functional differentiation analogous to (C14), in terms 
of mixed products of a M and fields 



(-0^:C(*i.-.«n) = jj^WpIK] 



= ( a% ( Xl ) . . . a% (x k ) A%+\ (x k+1 ) . . . A^ (x n ) >£> , (E2) 

K=0 



SWp - , , 

dW P 



6KW 

where the superscript (c) indicates the 'connected parts'. It follows then that 
where, for example, 

g {1 \ a {x) = {Al{x) )f + (al{x))f 
G (2 i?(x,y) = {A»(x)A b v (y))¥ + ( a a M (x)al{y) 
g^(x,y,z) = (A^x)A b v (y)Al(z))P + ( a a M (x)a b „(y)a%(z) , (E4) 

and similarly expressions higher order Green functions which involve 4, 5 ... space-time points. 



Wp of (El) involves the sources JC that do not have any immediate physical interpretation, it is more convenient 
to work with the corresponding effective action Tp, the generating functional for the proper vertex functions, which 
determines the equations of motion for the physically relevant Green functions. The effective action Tp is defined 
as the multiple Legendre transform, and is obtained by eliminating the source variables K in favor of the connected 
Green functions Q: 

Tp\g] = Wp[K] - /C«off« - l£(2)o(s( 2 ) 

_ 1^(3) ^(3) + 3 0(2)0(1) + 0(1)0(1)0(D) _ . . . . (E5) 

So far no approximations have been made. The variation of Tp with respect to the Green functions Q 1 ^ would yield 
an infinite set of coupled equations, the analogue of the BBGKY hierarchy |14[|: 



ST f 



5T 



£ = -i£( 2 )-i£(3)oaW 



gw 2 2 

^ = -^ (3) --' etc " ™ 

At this point approximation 1 and 2 of Section II A are invoked. It is assumed that the initial state is a (dilute) 
ensemble of hard gluons of very small spatial extent <C A, corresponding to transverse momenta k\ 3> /i 2 . By 
definition of A, or [i, the short-range character of these quantum fluctuations implies that the expectation value (a M ) 
vanishes at all times. However, the long-range correlations of the eventually populated soft modes with very small 
momenta k\ <C /x 2 may lead to a collective mean field with non- vanishing {A). Accordingly, the following condition 
is imposed on the expectation values of the fields: 
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= fort < t 
> for t > t 



( a a (x) ) = for alH 



(E7) 



Furthermore, the quantum fluctuations of the soft field are ignored, assuming any multi-point correlations of soft 
fields to be small, 



( A;\ ( Xl ) . . . AH (x n ) ) « < A% {x x ) > . . . ( A% (x n ) ) for all n > 2 



(E8) 



i.e., A^ is treated as a non-propagating and non- fluctuating, classical field. Hence, the set of Green functions (E4) 
reduces to 



G (1 V(x) = {A*Jx)) 



3» gW*{x,y) = {al{x)al{y))f 



(E9) 



These relations define the soft, classical mean field A, and the hard quantum propagators A. 

Now the hierarchy is truncated for n > 3. However, to perform this truncation properly, one must eliminate all the 
, etc., as dynamical variables by introducing p3] 



(E10) 



where 0(™- 1 for all n > 3 are functionals of of the 1- and 2-point functions alone, and are determined by the implicit 
equations 



;(") - 



0(1) 0(2) 



(2)1 



ST P [gw,g 
sgw [0W, 0( 2 )] 



= o 



for all n > 3 . 



(Ell) 



GfnJ = iA M „, so that Tp becomes a functional of only the soft mean field A^ and the hard propagator A^„: 



From (E3) and (E4) one sees that then the infinite set of Green functions reduces to involve only — A^ and 



rp [g] 



A, A 



= W P 



JC {1) o A - ~ /C (2) o f iA + A A 



(E12) 



The equations of m o tion for t he m ean field A and for the hard propagator A in the presence of external sources, 
follow now from (|e|), (|e|) and ( |eT2| ): 



<5r ? 



<5r P 



/CW^Cx) - / d*yKW^(x,y)A v '\y) 



5Af v (x,y) 2z 



£(2) Wi-fj;^) 



(E13) 
(E14) 



APPENDIX F: ANALYTIC PROPERTIES OF THE FREE-FIELD PROPAGATORS 

The components of the free-field propagator A " b are defined as in ( B19| ), i.e., 

Aft*(*,tf) = -i{Tal{x)a b M ) Aj^(x,y) = -i(T o«(x)o*(y) ) 

A^ b (x, y) = -i< a»a* (y) ) A<$(x, y) = -i( 4(vK(v) ) . (Fl) 

For free fields, one may write 

A $(x,y) = S ab d^(d x ) A (x,y) (A = A F , A F , A>, A<) , (F2) 



where d /iu (d x ) is defined by (j25|), and the functions Aq on the right side are the scalar parts of the propagators. The 
F,F,>,< components of the latter obey the following free-field equations with different boundary conditions: 
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-,2 
a 

O x 




^-2 


-\-o {x,y) 


-+2 
O x 


A^(r v) 

i x i y) 


— ^-2 

— A F (t h\ £> — 

— i x i y) Oy — 




->2 

d x 




= A>(x,y) dl = 







A* fay) 


^-2 

= A o( x >v) dy = 


0, 



and the identities 



(F3) 



&o(x,y) = 0(xo,yo) &o(x,y) + 0(y o ,x o ) A$(x,y) 

A F (x,y)=6(x ,y )A<(x,y) + 6(y , x ) A>(x, y) . (F4) 



Because of the relations (|F4|), the set of equations (F3) can be solved by only two independent functions, namely, (i) 
a purely imaginary and odd function i A~ , and (ii) a purely real and even function A + : 

iA-(x,y) = S ab d IM ,(d x )([a fl (x),ai(y)]) = i (A> - A<) (x,y) 
A+(x,y) = S ab M^MWO*) >«*(»)}) = i{^+A<)(x,y) (F5) 



From (F3) it follows that these functions obey 



dx A-(x, tf ) = A-(x,y) d y = A-(x,y) = -A~(y,x) 

dl A + (x,y)=A+(x,y) dl = A + (x,y) = +A+(y,x) , (F6) 



with the general solutions 

A-[x,y) = -i J -0^ e - ik <*-^ 2nS(k 2 ) ( 9l (k) g 2 (-k)) 

A+{x > y) = Sw? e ~ iHx ~ v) 27n5(fc2) + 52( ~ fc)) ' (F7) 

where the functions 

9l (k) = 8(k Q ) + h(k) g 2 (-k) = 6(-ko) + /a(-fc) (F8) 

contain the positive and negative frequency modes, respectively. Here 0(±fco) is the vacuum contribution, while 
fi t2 (±k) are the additional contributions from a medium. 

From ([F^)-(|F7|), one can now infer immediately the analytic properties of f\, f 2 , corresponding to those of <?i, g 2 . 

1. First, one observes, because A~ is purely imaginary and A + is purely real, that it must hold 

h(k)J 2 (k) = real. 

2. Second, because the commutator of free fields, i.e. the imaginary function A - , must be independent of the state 

of the medium, 

A-(x,y)= -i J ^ e -iH*-v) 2n6(k 2 ) (6(k ) - 6(-k )) , 

it follows that 

h{k) = M-k) = f(k) . 

3. Finally, because the anticommutator, i.e. the real function A + must satisfy 

d 4 xd 4 y h*(x) A+(x,y) h(y) > 



49 



for any smooth, but in general complex-valued function h, it follows that 

f(k) = fdk a f(k a ,k) > for all k , 



and so the 'on-shell' function f(k) is positive definite may indeed be identified with the positive definite phase- 
space density dN/d 3 k. 

The free-field solutions of A F , A F , A > , A< can now easily reconstructed using the following identities implied by 
© and (H): 



2A F (x,y) = -iA+(x,y) + (29(x° ,y°) - l) A~ (x,y) 
2A T (x,y) = -iA+(x,y) + (26(y°, x°) - l) A~(x, y) 



(F9) 



2A>(x,y) = -iA+{x,y) + A~(x,y) 

2A<(x,y) = -iA+(x,y) - A~{x,y) , (F10) 



from which, upon Fourier transformation, one obtains 

d 4 k _.-!./„_..•, / +1 



A F ix,y) = j— i e-^ ^ - i2^)f(k) 
A > Kj/) = J ^L e -*(-iO ( - i2 ^(fc 2 ) ((9(*b) +/(*))) 

A < (x,y) = | ^* c-*^-*) ( - z2^(fc 2 ) (0(-fc o ) + /(fc))) . (Fll) 

Finally, it is straight forward to infer the corresponding free-field forms of the retarded, advanced and correlation 
functions. 



A^(x,y) = +e(x ,y )A-(x,y) = (A F - A<) (x, y) = J ^e"*^) 
A^( : r,y) = -%o^o)A-(a ; ,y) = (A F - A>) (x,y) = | ^ e ~*(^) ( ' 



2ie 



A^fo^-iA+foy) - (A> + A<) (x, y) 

d 4 fc 
(27) 



-^ e -*(*-») ( - i 2 7 r<5(fc 2 ) (l + 2/(fc))) . (F12) 



50 



REFERENCES 



A. Billoire, G. Lazarides, and Q. Shan, Phys. Lett. B103, 450 (1 981); T. S. Biro and B. Miiller, Nucl. Phys. A561, 477 



(1993); G. Alexanian and V. P. Nair, preprint CCNY-HEP 2/95 (hep-ph/9504256) 



M. Bonini, M. D'Attana sio, and G. March esini, Nucl. Phys B444, 602 (1995); M. Bonini, G. Marchesini, and M. Simionato, 



preprint UPRF-96-464 (bep-ph/9604114) 



M. Baker, J. S. Ball, and F. Zachariasen, Nucl. Phys. B186, 531 (1981); M. Baker and J. S. Ball, Nucl. Phys. B229, 456 
(1983). 

See e.g., G. Marchesini, Nucl. Phys. 445, 49 (1 995); J. Kwiecins ki, Jour. Phys. G22, 685 (1996). 



R. Rajugopalan, preprint DOE-ER-40561-261, ( |trep-ph/9605327| ) 



M. Gyulassy and X.-N. Wang, Nucl. Phys. B420, 583 (1994); R. Baier, Yu. L. Dokshitzer, S. Peigne, and D. Schiff, Phys. 
Lett. B345, 277 (1995); E. Levin, in the proceedings of the RHIC Summer Study 96, edited by D. Kahana, Y. Pang and 



S. Kahana, BNL-publication (1996) (tiep-ph/9608304) 



J.-P. Blaizot and E. Iancu, in Quark-Gluon Plasma 2, ed. R. Hwa (World Scientific, Singapore 1996); J. -P. Blaizot and E. 
Iancu, to appear in Phys. Rep. 

R. D. Pisarski, Phys. Rev. Lett. 63, 1129 (1989); J. Frenkel and J. C. Taylor, Nucl. Phys. B334, 199 (1990); E. Braaten 
and R. D. Pisarski, Nucl. Phys. B337, 569 (1990); R. Kobes, G. Kunstatter, and A. Rebhahn, Nucl. Phys. B355, 1 (1991). 
A. Weldon, Phys. Rev. D42, 2384 (1990); E. Braaten and R. D. Pisarski, Phys. Rev. D42, 2156 (1990); ibid. D46, 1829 
(1992); P. Bozek and M. Ploszajczak, Phys. Rev. D41, 634 (1990); E. Braaten and M. Thoma, Phys. Rev. D44, 1298 
(1991). 

J.-P. Blaizot and E. Iancu, Nucl. P hys. B417, 608 (1993); ibid. B421, 565 (1994); ibid. B459, 588 (1996); preprint 



SACLAY-T96-085 ( [hep-ph/9607303; ), Saclay 1996. 
K. Geiger, Phys. Rev. D54, 949 (1996). 

Yu. L. Dokshitzer, V. A. Khoze, A. H. Mueller, and S. I. Troyan, Basics of Perturbative QCD, Advanced Series on Directions 
in High Energy Physics, Vol. 5 (Editions Frontieres, Gi f-sur-Yvette Ced ex, France, 1991) 



C. Greiner and B. Miiller, preprint DUKE-TH-96-99 ( |hep-ph/9605048| ) D. Bodeker, L. McLerran and A. Smilga, Phys. 

Rev. D52, 4675 (1995); F. Lombardo and F. D. Mazzitelli, Phys. Rev. D53, 2001 (1996). 

For a review see e.g.: H.-Th. Elze and U. Heinz, CERN Report TH-5325/89, 1989 (unpublished). 

J. Schwinger, J. Math. Phys. 2, 407 (1961). 

L. V. Keldysh, Sov. Phys. JETP 20, 1018 (1965). 

L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, New York 1962). 

E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics, Pergamon Press 1981. 

K. C. Chou, Z.-B. Su, B.-L. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 

E. Calzetta and B. L. Hu, Phys. Rev. D37, 2887 (1988). 

R. D. Jordan, Phys. Rev. D33, 44 (1986). 

J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

S. Mrowczyhski and U. Heinz, Ann. Phys.229, 1 (1994). 

L. McLerran and R. Venugopalan, Phys. Rev. D49, 2233 (1994); ibid. 3352; A. Ayala, J. Jalilian-Marian, L. McLerran 
and R. Venugopalan, Phys. Rev. D53 , 458 (1996); J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert, preprint 



TPI-MINN-96-1429 (tiep-ph/9606337) 



A. Makhlin, Phys. Rev. C51, 3454 (1995); ibid. C52 (1995). 

D. Boyanovski, H. J. Vega, R. Holman, D. S. Lee, and A. Singh, Phys. Rev. D51, 4419 (1994); D. Boyanovski, M. D. 
Attanasio, H. J. Vega, and R. Holman, Phys. Rev. D52, 6805 (1995). 
G. Leibbrandt, Rev. Mod. Phys. 59, 1067 (1987), and references therein. 

A. Bassetto, G. Nardelli and R. Soldati, Yang-Mills theories in algebraic non-covariant gauges (World Scientific, Singapore 
1991). 

B. S. deWitt, Phys. Rev. 162, 1195 (1067); L. F. Abbott, Nucl. Phys. B185, 189 (1981); D. G. Boulware, Phys. Rev. 
D23, 389 (1981); H.-T. Elze, Z. Phys. C47, 647 (1990) 

J. M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. D10, 2428 (1974). 

E. P. Wigner, Phys. Rev. 40, 749 (1932); M. Hillery, R. F. O'Connell, M. O. Scully, and E. P. Wigner, Phys. Rep. 106, 
121 (1984). 

[32] It should be remarked that the Wigner function can be transformed from the here employed class of non-covariant gauges 
to any other gauge by utilizing the general form of (|6^) with 

Gf v (r,s) = U aa '(r,r + ^j 0^' + |, r - |) U W (r - |, rj , 



51 



and U(x, y) is the 'parallel transporter' or 'link operator' (t],[16| that guarantees that the Wigner function transforms 
covariantly under local gauge transformations, when moving on an oriented straight line := x + £(y — x), < £ < 1 
from space-time point x to point y: 

ig J dz»A^(z)j = S ab ( 1 + P ^I1 fo*' J . 

The symbol P denotes the path-ordering of the color matrices in the exponential and A = T a A°^ where A stands for A 
or a. Choosing the gauge vector such that it is parallel to the direction of propagation, n \\ z, the gauge constraint 

n ■ A = implies z ■ A = 0, and as a consequence one has then U ab (x, y) " = ° S ab , so that 

y»v(r,k) = dse g M „ lr+ -,r - -I . 



This property is another very advantagous feature of the non-covariant gauges: the link operator becomes a unit matrix in 
color space, and so the complicated non-linear dependence (in the gauge fields) of the Wigner functions reduces to a form 
identical to those in classical field theory. This is of course not possible in a general gauge, in particular, in covariant gauges 
one has to account for the infinite number of non-linear terms in the series representation of the path-ordered exponential. 



See e.g.: J. D. Jackson, Classical Electrodynamics (Wiley, N.Y. 1975). 

Yu. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977); G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977). 
Yu. L. Dokshitzer, D. I. Dyakonov, and S. I. Troyan, Phys. Rep. 58, 269 (1980). 
A. Bassetto, M. Ciafaloni and G. Marchesini, Phys. Rep. 100, 203 (1983). 
K. Geiger, Phys. Rev. D51, 3669 (1995). 

A. Basetto, M. Ciafaloni, and G. Marchesini, Nucl. Phys. B163, 477 (1980). 
W. Furmanski, R. Petronzio, and S. Pokorski, Nucl. Phys. B155, 253 (1979). 
A. Basetto, M. Ciafaloni, G. Marchesini, and A. H. Mueller, Nucl. Phys. B207, 189 (1982). 
A. H. Mueller, Nucl. Phys. B228, 351 (1983). 

G. Marchesini, L. Trentadue, and G. Ve neziano, Nucl. Ph ys. B181, 355 (1981). 
Y. V. Kovchegov, CUNY-preprint 1996 jhep-ph/9605446| ). 
A. H. Mueller, Nucl. Phys. B307, 34 (1988). 

L. D. Fadeev and V. N. Popov, Phys. Lett. B25, 29 (1967); L. D. Fadeev and A. A. Slavnov, Gauge Fields (Benjamin 
Cummings 1980). 
[46] H.-T. Elze and U. Heinz, Phys. Rep. 183, 81 (1989). 



52 



